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Abstract 

The transverse spatial effects observed in photon pairs produced by paramet- 
ric down-conversion provide a robust and fertile testing ground for studies 
of quantum mechanics, non-classical states of light, correlated imaging and 
quantum information. Over the last 20 years there has been much progress in 
this area, ranging from technical advances and applications such as quantum 
imaging to investigations of fundamental aspects of quantum physics such 
as complementarity relations. Bell's inequality violation and entanglement. 
The field has grown immensely: a quick search shows that there are hun- 
dreds of papers published in this field, some with hundreds of citations. The 
objective of this article is to review the building blocks and major theoretical 
and experimental advances in the field, along with some possible technical 
applications and connections to other research areas. 
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1. Introduction 

Correlations between optical fields produced in spontaneous parametric 
down-conversion (SPDC) were first observed in 1970 by Burnham and Wein- 
berg [l[. In this early experiment, they investigated the temporal and spatial 
correlation between pairs of photons. From the theoretical point of view, it 
was Zeldovich, Krindach and Klyshko who first studied the statistical proper- 
ties of the light produced in this process, as early as 1969 ^,i3||. The temporal 
correlations were studied and utilized in the pioneering series of work of Prof. 
Leonard Mandel and his group 0, Is], @, 0, Isl , giving rise to the use of the term 
twin photons, since the temporal correlation indicates that pairs of photons 
are born simultaneously. These correlations are of fundamental importance 
in all experiments with twin photons. They provide the time correlation 
which allows for the post-selection of photon pairs born from the same pump 
photon, which may present correlations in other degrees of freedom, such 
as the transverse position and momentum. Both the time and spatial cor- 
relations have been proven to be non-classical. The simple observation of 
a coincidence rate above the rate of accidental coincidence counts already 
leads to the violation of a classical inequality . In the early Burnham and 
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Weinberg experiment spatial correlations were already observed. They 
demonstrated that the transverse profile of the coincidence distribution in 
the far field is narrower than the single photon intensity distribution. How- 
ever, more comprehensive investigations were made only in the 1990s. For 
instance, Grayson and Barbosa 10| investigated spatial correlations for use 
in induced coherence without induced emission 11. Joobeur et al. 12 stud- 



ied general spatial and temporal coherence properties, and Souto Ribeiro et 



al. [13(1 started a long series of experiments using double-slit apertures and 



coincidence detection, initially dedicated to the study of spatial coherence 
properties of twin photons. Strekalov et al. 1J| performed a novel double-slit 
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ini- 



experiment to observe ghost interference fringes and Pittman et al. 
tiated a series of experiments on what was called Quantum Imaging, with 
the two-photon field from SPDC Monken et al. 16| demonstrated that the 
spatial correlations could be controlled through the shaping of the angular 
spectrum of the pump beam. This control allowed the use of the spatial 
correlations in a series of applications, ranging from fundamental aspects of 
Quantum Mechanics such as the measurement of the photonic de Broglie 



wavelength of a two-photon wave packet |17l | and the observation of the spa- 



tial anti-bunching of photons [18|,ll9|], to applications to quantum information 
such as, for instance, the production of spatial qudits [201 . Howell et al. |2i| 
and D'Angelo et al. 22| formally demonstrated entanglement between the 



spatial properties of the twin photons. 

In this review, we will discuss the spatial correlations between twin pho- 
tons produced in parametric down-conversion, starting with the basic con- 
cepts and basic theory for the description of the down- conversion process 
and temporal simultaneity. We will introduce some basic concepts related 
to the classical spatial coherence of light fields and propagation of parax- 
ial light beams. The theoretical framework for the description of quantum 
spatial correlations will be introduced and the remainder of the text will 
be dedicated to the description and discussion of several key experiments 
and applications: double-slit experiments, quantum images, demonstration 
of non-classicality, spatial entanglement, transverse modes and applications 
to quantum information. 



2. Fundamentals of parametric down-conversion 

In SPDC, a nonlinear birefringent crystal is typically pumped by an in- 
tense pump laser beam, producing low intensity signal and idler fields, as 
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Idler 




Figure 1: Pumping a non- linear crystal with an intense laser beam produces low inten- 
sity signal and idler beams by the non-linear process of spontaneous parametric down- 
conversion. 



illustrated in Fig. [TJ The pump laser has frequency Up and is typically in 
the UV or violet region, while the down-converted signal and idler fields have 
frequencies Ug and Ui that are usually in the visible or near infra-red region 
of the spectrum. There have been many theoretical treatments of SPDC, 
beginning with the early work of Klyshko [233 later by Hong and Mandel 
A detailed account of these calculations is out of the scope of this review, 
but can be found in the Ph.D. thesis of L. Wang, 2^ and the books by D. N. 



Klyshko (25|, Mandel and Wolfl26|Land Ou|27 
treatment presented in Refs. |5|, |2; 



In this section, we follow the 



to introduce the main features related 
to the temporal correlations between twin photons produced in SPDC. Even 
though our main concern is the spatial correlation, their time correlation is 
essential in all experiments. 

2.1. Classical and quantum parametric interaction 

The quantum theory of parametric down-conversion can be derived from 
the classical description of the nonlinear interaction, followed by the quan- 
tization of the electromagnetic field. Up to second order, the components 
of the electric polarization of a nonlinear and non centrosymmetric optical 
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medium when an electric optical field E propagates through it is [28|, |29 

oo 

P,(r,t) = eoJdt'x^\t')E^ir,t~t')+ (1) 



oo oo 

I dt' I df'x'^lit', t")E,{r, t - t')E,{r, t - t"), 







where Ej{r, t) is the j component of the electric field vector and x^^^ and x''^'* 
are the first and second order susceptibility tensors, respectively. Summation 
over all combinations of the tensor components and electric field components 
is assumed. For light sources with low electric field strengths compared with 



the characteristic atomic electric field strength Eat ~ 6 x 10^^ V/m 30|], only 
the linear term in Eq. ([1]) is significant. Increasing the field intensity by using 
laser sources, for example, increases the nonlinear terms in the polarization. 
We will focus on the second order non-linear interaction. 

In order to quantize the field, we begin with the electromagnetic field 
Hamiltonian in the dielectric medium of volume V / citemandel95 

n{t) = \j dv [D(r, t) ■ E(r, t) + B(r, t) ■ H(r, t)] , (2) 

V 

where D is the displacement vector, B is the magnetic induction and H is 
the magnetic field. If we use the definition D(r, t) = eoE(r, t) + P(r, t) in 
Eq. 02]), and use Eq. ([T]), we can rewrite the Hamiltonian as 

'H(t)=Ho{t)+Hi{t), (3) 

where l-Loit) contains the interaction of the electric field and the first order 
linear component of the electric polarization. The "perturbation" 'Hi{t) is 
the nonlinear interaction Hamiltonian, given by 

•H/(t)=^ J t^rE(r,t)-P,i(r,t) 

V 

oo oo 

= \j dv j dt' j df'x'^lit', t")E,{r, t)Ej{r, t - t')E,{r, t - t"), (4) 



V 
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where Pni(i", t) is the nonhnear component of the electric polarization, and 
summation on repeated indices is understood. In order to avoid unnecessary 
difficulties in the quantization of the field, it is considered that there are no 
electromagnetic boundaries between the medium and the air. In all situa- 
tions discussed here, the effects related to refraction and birefringence can be 
included classically, after the field is quantized. A more rigorous approach is 
not yet available and is still subject of research. 

From this point on we concern ourselves only with the nonlinear portion 
of the Hamiltonian. It is also convenient to suppose that just after the crystal 
there are two interference filters limiting the frequency spectrum of the down- 
converted fields, and that the pump beam has a narrow frequency spectrum. 
The two-photon quantum state will then include the functions describing the 
signal and idler interference filters. As a part of the quantization procedure, 
we expand the classical optical electric field in terms of plane waves: 

E(r,t) = E+(r,t) + E-(r,t), (5) 

with 

E+(r, t) = ^ y2 ek,a£k,.o;k,aG'(a;) exp [i(k • r - ut)] = [E" (r, t)] * , (6) 

where 

£k,a = \/fioj{k, (7)/2eon2(k,(7), (7) 

eo is the free space permittivity, G{uj) is the transmission function of the 
interference filter, V is the quantization volume, k is the wave vector, u is 
the two-dimensional polarization vector, u is the frequency, and Q;k,cr is the 
mode amplitude. The index a is summed over orthogonal components of a 
two-dimensional polarization vector and k is summed over all possible wave 
vectors. 

We adopt the usual method of quantization of the electric field, letting 
ak,o- — > ak,crj where ak,o- is the photon annihilation operator. Then, the 
electric field amplitude becomes a field operator, given by 

E+(r, t) — > E+(r, t) = ek,a£k,aak,aG(cc;) exp [i{k ■ r - ut)\ 

= [E-(r,^)]^ (8) 
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Substituting expression into the classical Hamiltonian (jl]) we have a quan- 
tum Hamiltonian operator 

^ 2y3/2 X] 3t,,as 9t,,a, 9kp,apai^^^^ai^^a,^kp,ap 

l^CTi'^S kjjCTj kp,f7p 

X exp [i{us + uJi - Up)t] Xijfc(ek„aJ* (ek„aj*(ekp,ajfe 

X j exp [-iiks + ki- kp) ■ r] rfr + H.C., (9) 
y 

where we have distinguished the three fields as pump (p), signal (s) and idler 



*^^--H/^^^^^GMk,..,)l. (10) 

V is the interaction (crystal) volume, and H.C. stands for Hermitian Conju- 
gate. Here, ^(kj, aj) is the linear refractive index of the (anisotropic) crystal. 
We have also eliminated all terms which do not conserve energy and have 
defined 

Xijk = x!J(Wp = Us + UJi)+ x'^^li^i = ^s+ Wp) + xJJ(Ws =UJi+Up), (11) 

with 



oo oo 

xgi(a; = uj' + uj") = jdt' j df'x^lit', t") exp [-^{uJ't' + J't")]] . (12) 





To find the quantum state at a given time t, we assume that the nonlinear 
interaction is turned on at time to = when the system is in the initial state 
|'?/'(0)). The state at time t is given by the time evolution of some initial state 
at to = 0: 

m)) = m \m) (13) 

where 

U(t) = exp j rfrH,(r)j , (14) 
is the time evolution operator. 
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If the pump field is sufficiently weak, such that the interaction time is 
small compared to the average time between down-conversions, then we can 
expand Eq. f|T^ in power series and keep only the two-photon term: 

U(t) = l+ j (irH,(r)j +■■■. (15) 

The integral can be expressed in terms of Eq. ([9]) as: 

ks,crs kjjfJ^ kp,f7p 

X exp [i{ujs + uji - ujp)t/2\ Xijfc(ek„aJ • (ek,,<7j*(ekp,ajfc 
X tsinc[(a;s + oji — Up)t/2] j dr exp [— i(ks + kj — kp) ■ r] 

V 

+H.C. (16) 

Integration in r leads to a sine function involving the wave vectors, which 
provides the conservation of momentum condition: 

ks.fJs kijfjj kp,crp 

X Xijfc(ek„aJi (ek,,aJK^kp,aJfc sinc[(w, + Ui - Up)t/2\ 
X exp [i{uJs + ujp)t/2\ J]^ sine [(k^ + kj - kp)m/m/2] 

m 

X exp [-z(k, + k, - kp) J,/2] + H.C., (17) 

where V = Ix ^ ly x Iz and Im is the dimension of the nonlinear medium in 
direction m {m = x, y, z). 

The quantum state at time t can be finally obtained using Eq. ( fT3l) . 
considering the initial state as the vacuum state, and using the interaction 
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Hamiltonian in the form of Eq. f ll7l) : 



|^(t)) = \vac) + ^^^^3^, sL-s 9L<r,. ^7k„.,^^p(kp, ap) 

ka,f7s ki,0'i kpjCTp 

X Xiifc(ek„<xJ* (ek„aj*(ekp,ajfcsinc [(w^ + - Wp)t/2] 
X exp [i{us + Up)t/2] J]^ sine [(k^ + kj - kp)„/„/2] 

m 

X exp [-z(ks + ki - kp)Jj2] \ks, as) |ki, aj) , (18) 

where |ks, cTs) and |kj,(Tj) are single photon Fock states in modes (kg, as) 
(signal) and (kj,(Tj) (idler) respectively, and fp(kp,crp) is a classical ampli- 
tude corresponding to the plane wave component (kp, Up) of the pump beam. 
Quantum effects in the pump field are neglected in this treatment, since the 
pump depletion is much smaller than its average intensity and its quantum 
state is effectively constant. The annihilation operator for the pump modes 
was replaced in Eq. f|T8l) by a classical amplitude. 

In order to simplify expression f[T8l) , it is convenient to make the following 
approximat ions : 

(a) The interaction time is long enough, so that the term sine [us + Ui = 
cOp)t/2 is significant only when ojs+coi = Up. This assumption can be justified 
by the use of a moderate power pump laser so that the time interval between 
two down-conversions is large compared to the detection resolving time. 

(b) The frequency spread of the detectable down-converted fields is small 
compared to the central frequencies, so that the dispersion of the refractive 
indices around the central frequencies Uj is small and a linear approxima- 
tion can be used. This assumption is justified by the use of narrow-band 
interference filters in front of the detectors. 

~ (2) 

(c) The terms fi'kj.o-j and xlji slowly- varying functions of k^, so that 
they may be taken as constants in the intervals considered for kj. 

(d) The pump beam propagates along the z axis and the crystal is large 
enough in the x and y directions to contain the whole pump beam transverse 
profile. In this case, Ix and ly can be extended to infinity and the last term 
in the third line of expression f|T8l) is proportional to 

^(q. + q* - %) sine [{ksz + ki, - kp^)L/2], 

where = [kj^, kjy) is the transverse {xy) component of kj and L = is 
the crystal thickness. 
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(e) The quantization volume is large enough to justify the replacement 
of summations over k by integrals. 

(/) The pump beam contains only extraordinary polarization. It is im- 
plicit in this assumption that we are dealing with negative birefringent crys- 
tals. 

Under the above assumptions, Eq. f lTB]) is written as 



\vac 



) + ^ j dus j dui j dq^ j dq, (q„ q^, w„ Wj) 

X |q,,cUs,o-^) \(l„Ui,a^) , (19) 

where \q^j,cuj,aj') represents a one-photon state in the mode defined by the 
transverse component q^- of the wave vector, by the frequency Uj and by the 
polarization aj. The amplitude $ is now reduced to 

^asa, ~ CcTsa, Gsiujs)Gi{ui) f (q^ + q^, + Wi) sluc [{ksz + ki;,-kp^)L/2], (20) 

where C^^^. is a coupling constant which depends on the nonlinear suscep- 
tibility tensor, and G{ujj) is the spectral function defined by the narrow 
bandwidth filters placed in front of the detectors. If the anisotropy of the 
medium is neglected (which is not always convenient), the longitudinal wave 
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vector mismatch ksz + — kpz is written as y^|k7p"-^]q7P + y^lkiP — |qj 

2.2. The coincidence count rate 

The state given by Eq. f llSp can be simplified, considering approxima- 
tions that are appropriate to given specific experimental situations. In this 
section, let us consider that signal and idler fields are detected through small 
apertures at fixed positions far enough from the source, so that only one spa- 
tial mode is selected by each detector. In this case, when narrow bandwidth 
filters are used, the quantum state takes a very simple form: 

\i){t)) = Ci \vac) + ^2 j dus j dui G s{i^s)Gi{ui)vp{uJi + Us) \uJs) , (21) 

where Ci and C2 are normalization constants. 

The electric field operator can also be simplified to 

E+(t + r) = C j da;a(a;)exp[-ia;(t + r)], (22) 
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where C is a constant. 

The state in Eq. ( ETj) and field operator in Eq. (12^ were obtained under 
several approximations. However, they adequately describe the main features 
of time correlations between twin photons. The time correlation is always 
assumed and used in all twin-photon experiments with coincidence detection, 
even when the spatial correlations are the focus of the investigation. There- 
fore, it is important and interesting to calculate the coincidence count rate 
in this simplified point of view. 

The coincidence count rate is given by the fourth order (in the fields) 
correlation function: 

Rcit + Tut + Ts) = (E;(t + r,)Er(t + r,;)E+(t + ri)E+(t 



r. 



)) 



(23) 



Using Eqs. ([2TD and ([22]), we obtain: 



Rc(t + Ti,t + Ts 



j dui Jdu2 a(wi)a(a;2) exp {-i[ui(t + ti) + U2{t + r^)]} 



X dui dusGi{uji)Gs{us)vp{ui + ujs) \ujs) 



(24) 



where 77 is a constant that depends on the square root of the pump beam 
power, the efficiency of the detectors and of the parametric down conver- 
sion process. After the action of the operators on the quantum state and 
integration over uji and a;2, we obtain: 



V 



doji dujsGi{uji)Gs{ujs)vp{uji + ujs) 



X exp + uJs)t] exp [-i{uiTi + UsTs)] 



■ (25) 



In the limit where the pump spectrum can be approximated by a delta 
function Vp{u) — )■ 6{up — u), the time correlation is: 

2 



Rc{Ti,Ts) 



V 



duj Gi{uj)Gs{ujp — uj) exp [iuj{Ts — Tj)] 



(26) 



where F{t) is the convolution of the Fourier transforms of the filter functions 
Gi{(jj) and Gs{ujp — u). Usually, the width of F{t) is on the order of fem- 
toseconds. This means that detections are simultaneous within the window 
of a few femtoseconds. 
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3. Fundamentals of spatial correlations 



The investigation of tlie transverse spatial properties of light has accom- 
panied the study of the nature of light itself. The development of the theory 
of diffraction and all resulting applications represent an impressive example 
of how finely a physical theory can describe nature. Despite the incredible 
success of classical optics, it has been found in the last few decades that 
amendments should be made in order to accurately describe the quantum 
effects which arise when dealing with some special light sources. In this 
respect, the transverse correlations between the twin photons produced in 
the parametric down-conversion process have largely contributed to this end. 
The quantum optics theory of spatial effects which has been developed so 
far is strongly based on classical diffraction theory, and has been shown to 
describe a large array of quantum phenomena in optics. In the following, we 
will briefly discuss some important topics in classical diffraction theory which 
will help in the development of the quantum theory of spatial correlations. 

3.1. Transverse coherence and partial coherence: classical optics 

In this section we will review some fundamental aspects of transverse 
coherence and partial coherence in classical optics. These concepts were es- 
sential in the development of quantum coherence theory and are still essential 
for the design and understanding of many quantum optics experiments. We 
introduce the discussion by analyzing the Young's double-slit experiment and 



its relationship with the van Cittert-Zernike (3l|] theorem, and the concept 
of coherence area. 

3.1.1. Double slit experiment: perfect coherence 

Fig. [2] shows a light source illuminating a double-slit aperture, and the 
resulting intensity pattern which is observed at the detection screen. When 
the light emitted by the source is perfectly coherent, interference fringes are 
observed with high visibility or contrast, as in Fig. [2^). Many of the laser 
light sources available today are almost perfectly coherent, but many other 
important light sources are incoherent or partially coherent. Incoherent light 
sources can also be used to perform interference experiments under certain 
conditions. All real light beams present partial coherence, depending on 
the geometrical properties and propagation of the beam. Thus, interference 
fringes can still be observed, but with reduced visibility, as illustrated in Fig. 



Wp)- The light intensity distribution after a double-slit is given by |32l. |31 
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a) Jk. b) J\f^ 




Figure 2: Double-slit experiment. The double-slit is illuminated by a perfect coherent 
light source in a) and by a partially coherent light source in b). 



I{x,y) = Io{x,y)[l + \-fi2{d)\cos(k- p + 6)], (27) 

where /o(a;, y) gives the single sht diffraction pattern, k is the wave vector, p 
gives the position in the transverse plane, 712 is the normalized mutual coher- 
ence of the light field in the plane of the slits, d is the distance between slits 
and (5 is a fixed phase factor. When the source is perfectly coherent, I712I = 1 
and the contrast of the interference fringes is maximum. The contrast can 
be quantified by the visibility V, defined as = (/max — /min)/(/max + Imm), 
where /max and /min are the intensities at the interference maximum and min- 
imum. When the coherence is partial, < I712I < 1, and can be calculated 
from geometrical and statistical properties of the source 31 . 

One important result in the calculation of 712 is known as the van Cittert- 
Zernike theorem. We will discuss this theorem in more detail. 

3.1.2. Partial coherence and the van Cittert-Zernike theorem 

Let us consider the situation sketched in Fig. |3l where a light source 
S is illuminating the observation screen O. We assume that the source is 
constituted by infinitesimal independent light emitters, such as a thermal 
source or parametric down-conversion, for example. The non-normalized 
mutual coherence function between points 1 and 2 on the screen is defined 
as: 

Tuid, t) = J2{Emi{t + r) EUt)) + E E + ^) W), (28) 

where the summations are performed over the emitting points of the source, 
d is the distance between points 1 and 2 at the screen, r is the time difference 
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Figure 3: Finite source S illuminating an observation screen O. 



due to different propagation distances from point m at the source to point 
l(or point 2) at the screen, E* is the complex conjugate of E, and the brackets 
indicate time average. 

If the field is quasi-monochromatic, the mutual coherence function can 
be written as 



where k = |k| and c is the speed of light. This approximation is valid when 
the time difference r is much smaller than the coherence time of the source. 
In other words, the phase relationship between Emiit + t) and Em2(t) is 
preserved. However, for spatially incoherent sources such as thermal sources 
and SPDC, there is no phase relationship between Emiit + r) and En2{t) if 
m ^ n and the second term in the left side of Eq. (I28l) averages to zero. 
Therefore, we end up with: 



It is convenient to take the continuum limit and suppose that each point 
of the source emits a spherical wave. Thus, Eq. (!30|) takes the form: 




(29) 




(30) 



m 




(31) 



s 
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R 



Figure 4: Finite source illuminating an observation screen. 



where /(ro) is the intensity distribution of the source, Ri = {tq — ri| and 
R2 = 1 1*0 — 1*2 1 are the distances between point tq in the source and points ri 
and r2, respectively, and it was assumed that Ri R2 = R ^ {Ri — R2)- 

It is interesting to write the difference Ri — R2 exphcitly in terms of 
the coordinates of the generic source point (xo,?/o) and the coordinates of 
the points 1 {xi,yi) and 2 {x2,y2) (see Fig. H]) and use again the relation 
i?i ~ -R2 = -R ^ {Ri — -^2)- We then obtain the result known as the van 
Cittert-Zernike theorem: 

Tuid, ^) = %jj dxodyo I{xo, yo)e''^'''°^'''\ (32) 
s 

where a is a constant phase and 

{xi - X2) [yi - 1/2) 
p= andg = . (33) 

The dependence on the distance d appears through p and q. The van Cittert- 
Zernike theorem shows that the mutual coherence is given by the Fourier 
transform of the intensity distribution of the light source. This result is valid 
if the emitting points in the source are independent and if the mutual coher- 
ence is considered in a plane far enough from the source, so that Fraunhofer 
diffraction regime can be assumed. The normalized degree of coherence 712 
is: 

e-//dxo%o/(xo,^o)e^^(P-°+''^») 
jj dxodyoI{xo,yo) 

The quantity |7i2('^)| gives the visibility of the interference fringes in a double 
slit experiment, Eq. ( l27l) . where we remember that d is the distance between 
slits. 



16 



Figure 5: Gaussian source illuminating an observation screen. 



3.1.3. The coherence area 

Let us illustrate the application of the van Cittert-Zernike theorem in 
the calculation of the coherence area of a light source. Referring to Fig. 
[5l we start by calculating the degree of coherence for a light source with 
a Gaussian intensity distribution. The integrals in Eq. (IMll are computed 
using the initial intensity distribution 



The coherence area is the area within which the value of I712I is appreciable 
and the field can be considered spatially coherent. From Eq. ( I36p . we can 
see that, on the one hand, the degree of coherence decreases with the wave 
number fc, source width a and distance d between points 1 and 2. On the 
other hand, the degree of coherence increases with the propagation distance 
R. As the coherence function is Gaussian, it never reaches zero and we can 
take the value 712(^^,0) = 1/e as a lower bound for the coherence. Then 
the circular coherence area Ac = 7r(|)^ can be deduced from Eq. (!36l) by 
imposing 712(6?, 0) = 1/e : 



In terms of a double slit experiment, this result indicates that there will be 
interference fringes with considerable visibility as long as the slits are placed 




(35) 



and we obtain 




(36) 
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within the area given by Eq. (1H7|) . 

The double-sht experiment provides a method for determining the coher- 
ence area of a hght field at some distance from the source. A double-sht 
aperture with variable slit separation d is placed at the plane of interest, the 
visibility of the resulting interference patterns is measured in the far-field. 
By plotting the visibilities of the interference patterns as a function of c?, a 
measurement of the transverse coherence length of the field in the one direc- 
tion is obtained. Repeating the procedure in the perpendicular direction in 
the transverse plane, we obtain a measurement of the coherence area. This 
method was used to measure the transverse coherence length of one of the 
down-converted fields produced by SPDC ■S^], and it was observed that the 
coherence length is the same as that predicted for a thermal source. The 
same method was also applied to the two-photon beam produced by SPDC, 
and it was shown that the transverse coherence length of the two-photon field 



is much larger than the one-photon field [3J] . A quantum multimode theory 



was necessary to explain the experimental results in this case. 

3.2. The Fresnel and paraxial approximations 

Another important and basic issue in spatial correlations concerns the 
Fresnel and paraxial approximations. Most of the quantum optics theory 
concerning SPDC was developed within these approximations and describes 
very well the vast majority of the experimental situations reported so far. 
In the following sections, we will derive the quantum state for the trans- 
verse degrees of freedom of the twin photons produced by parametric down- 
conversion, according to the Fresnel and paraxial approximations. 

Let us assume that a monochromatic light beam is propagating along the 
z direction in an isotropic medium. We can write the component as 

K=vi^^'^k(^i-^y (38) 

The above approximation, known as the Fresnel approximation, is obtained 
by simple Taylor expansion and is valid when ^ /e^. 

The Fresnel approximation is essentially a particular application of the 
more general paraxial approximation. In geometric optics, where light is rep- 
resented by rays, paraxial rays are those that lie at small angles to the optical 
axis of the optical system under consideration. If we were to draw rays from 
the origin to points k in /c-space satisfying the approximation in expression 
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they would be paraxial rays. In this respect the paraxial approxima- 
tion and the Fresnel approximation are essentially the same. Throughout 
this work we will refer to both as the paraxial approximation. A thorough 
account of the paraxial approximation can be found in most classical optics 



textbooks [35|,|36[. Here we will present only the bare essentials. The parax- 
ial approximation can be extended to wave optics if one considers only waves 
whose wavefront normals are paraxial rays. In wave optics, an optical wave 
E{r,t) satisfies the wave equation: 

V^E(r,t)-i^E(r,t) = 0. (39) 

If one considers that the optical wave is monochromatic with harmonic time 
dependence, so that E{r,t) = £{r) exp^—icot), where u is the angular fre- 
quency, one arrives at the Helmholtz equation: 

V^^(r) + k^S{r) = 0. (40) 

If we now consider only paraxial waves propagating near the z axis, we can 
write S{t) = U{t) ex]i{ikz), where W(r) is a slowly varying function of r such 
that £^(r) maintains a plane wave structure for distances within that of a 
wavelength. Using this form of ^^(r) in the Helmholtz equation, one arrives 
at the paraxial Helmholtz equation 

+ -^ + 2tk-^]Uir) = 0. (41) 



In obtaining fHT]) . we have used the fact that the term d'^U{r)/dz'^ is very 
small within distances of a wavelength: d'^U{r)/dz'^ <^ kdU{r)/dz. The 
paraxial Helmholtz equation admits several well known solutions, including 
the Hermite-Gaussian and Laguerre-Gaussian beams. Quantum fields pre- 
pared in these modes have been proposed for several applications and will 
be discussed in more detail later. It has been shown by several authors that 
the paraxial Helmholtz equation is analogous to the Schrodinger equation in 



quantum mechanics [37|, |38|, |39| . In reference [37[ an alternative derivation of 
equation fHTl) is provided which requires that the optical wave is only nearly 
monochromatic. 

3.3. The angular spectrum and its propagation 

The quantum theory which accounts for the transverse correlations of 
photons makes use of several techniques of Fourier Optics, in particular the 
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propagation of the angular spectrum. Fourier Optics [36[ provides a useful 
method of calculating the propagation of an optical field through a given 
optical system. Here we consider a monochromatic scalar field far from its 
source, satisfying the Helmholtz equation f HOj) . for which the field amplitude 
can be represented [36i] as 

W(p, z) = j^-^ j dq v{ci, z) exp [iq ■ p] , (42) 

where we assume that the field is propagating along the z direction and have 
defined p and q as the transverse components of r and k, respectively. The 
angular spectrum f (q, z) is the inverse Fourier transform of the electric field: 

One can also understand the angular spectrum by recognizing equation ( l42|l 
as an expansion of W(p, z^ in terms of plane waves exp(iq ■ p), in which the 
angular spectrum t>(q, z) acts as a weight function. 

As the field propagates, its angular spectrum changes accordingly. Com- 
bining Eqs. fHTj) and (H2|) . it is easy to show that if the angular spectrum of 
a field is known at z = 0, it propagates as 

f (q, 2;) = t>(q, 0) exp(iA;22;). (44) 

In order to be useful, must be expressed in terms of the transverse com- 
ponent q. For isotropic media, this can be done using Eq. (1551) . 

Dealing with anisotropic media, which is the case of the nonlinear crys- 
tals used for parametric down-conversion, the expressions for kz in terms of 
q can be quite complicated 32] • We will restrict our analysis to uniaxial 
media for two reasons: first, in the majority of the work on spatial corre- 
lation properties of two-photon states, the photon pairs are generated by 
spontaneous parametric down-conversion in uniaxial crystals. Second, the 
physics of down-conversion in biaxial crystals, though more involved, to our 
knowledge does not contain any essentially new effect. We will consider the 
geometry depicted in Fig. [6l where a slab of a uniaxial crystal of thickness 
L is cut to have its optic axis making an angle Q with its normal, which is 
oriented parallel to the z axis. In this medium, the wave vector surface has 
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Figure 6: Non-linear crystal and its optical axis, z axis represents the direction of propa- 
gation of light. 



two sheets, defined by 



ql + ql + kl uj^ 



^2 ^2 



(45) 



which allows us to immediately write 



and 



kz = \l { no-) - |qP ^ rio- - tt^Iq^, (46) 
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cos^ Q 


^ sin2^\ 


nl 


nl ) 




1\ . 




^ sm 




nl) 



+ ^ + 1 T- + 2 

ni \ ni ni 
'iJ_ 
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(47) 



where n^ and n^. are the ordinary and extraordinary refractive indices, respec- 
tively. Eq. ( H5|) applies to plane waves with ordinary polarization (orthogonal 
to the optical axis), whereas Eq. ( |47|) applies to plane waves with extraordi- 
nary polarization. In the latter case, Eq. ( H7|) can be solved for A;^, leading 
to 



2 

1 

y 



kz = aqx + y (v'^) - l^'^ql - 1^% 



+ (48) 
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where 



a 



/3 



7 



7] 



(nj, — nl) sin 6* cos 6^ 



n1 sin 6 + n1 cos^ 6 
n1 sin^ 6 + n1 cos^ 6 ' 
\/ nl sin^ 6^ + cos^ 9 



(49) 
(50) 
(51) 
(52) 



\/ n1 SIP? 9 + n1 cos^ 9 

The term a is responsible for the so-called walk-off. The terms /3 and 7 are 
close to 1 and cause a slight astigmatism in the beams propagating through 
the uniaxial medium. Their effect is marginal, and both and 7^ can be 
approximated by 1. Therefore, for extraordinary polarization, 

u c , ,0 



k^K.r]- + aq^ - - — |q| 
c irjuj 



(53) 



3.4- Quantum state of photon pairs: spatial degrees of freedom 

The spatial properties of the two-photon state depend strongly on the 



birefringence of the non-linear crystal and type of phase matching [40|, |41 



42I . I43I . |44| . Let us consider first the case of type I phase matching, where 



the pump field has extraordinary polarization and the down-converted fields 
have ordinary polarization (e — i- 00). In this case. 
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(54) 
(55) 
(56) 



where nj {j = s,i) stands for the ordinary refractive index. 

Using Eqs. (!54|) - (l56|) in Eq. (120|) and considering a monochromatic pump 
beam, we arrive at the two-photon detection amplitude for type I phase 
matching: 



[(Is, 



Coo Gs{ojs)Gi{ui) v{q^ + qi)S{us + Ui - 



X smc 



/ioo + k{qsx + Qix) + 



AK 



Si 



X fiy.Y>[-ilt{qsx + qix)], 



(57) 
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where K = rjpUJp/c is the pump beam wave number inside the crystal, 
\JuJs/uJu Ti = y/uj/us, jion = {u - r]p)KL/2r]p, and k = apL/2 is 



the transverse walk-off length^^ 43[. The ordinary refractive index n is cal- 
culated at the frequency Up/2. In collinear phase matching, = 0. Eq. 
( I57|) is accurate up to first order in {us — uJi)/2ujp. 

For each particular set of frequencies w^, Wj, the function $oo(qs,qj) can 
be regarded as the normalized angular spectrum of the two-photon field. Note 
that the angular spectrum of the pump beam is present in $00, meaning that 
it is transferred to the fourth-order spatial correlation properties of the two- 



photon state |16|. The amplitude ^00 is not a separable function of q^, and 
qj, that is to say, $oa 7^ Fs{.^.s)Fi{.^.^■ This non-separability is responsible for 
many of the nonlocal and non-classical effects observed with the state fl60l) . 

For simplicity, let us assume that the crystal is thin enough to allow 
us to ignore the effects of birefringence. If narrow-band interference filters 
selecting Ug = uji = ujp/2 are used, the two-photon state generated by SPDC 
in quasi-coUinear phase matching (/ioo ~ 0) can be written in a very simple 



form as |16l . |46 
where 



lV')spDC = '^v \vac) + C2 \ip)i , (58) 

■ j I dqjqi w(q, + qihiq^ - qj Iq^, 60) |q„ e^) , (59) 

7(q) = 2L /tc^ Ksmc{Lq'^ / AK) , and |qj,eo) represent Fock states in plane 
wave modes labeled by the transverse component q^- of the wave vector k^, 
and ordinary polarization vector Sq. If the nonlinear crystal is sufficiently 
thin, so that the width of the sine function in Eq. fl57p can be much greater 
than the width of the pump beam angular spectrum, the sine function can 
then be approximated by unity. This is known as the thin crystal approxi- 
mation. The quantum state f l59|l simplifies to 



1^) 



tc 



j j dqjqi v{q, + qj |q„ e^) [q^, e^) , (60) 



D 



D is a domain in the q space within which the thin-crystal approximation 
is valid, meaning that the two-photon angular spectrum $oo(q^, qj mimics 
the pump beam angular spectrum v{q) in the sum of transverse wave vectors 
q^ -|- qj. Since all the information about a monochromatic beam is contained 
in its angular spectrum, the transverse and longitudinal spatial properties 
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of the pump beam are transferred to the two-photon field, in this context, 
defining a correlation beam. 

In the general case, however, when the crystal length is not negligible, 
the sine function in Eq. fl57|) has to be considered in detail. Due to the 
presence of the linear term lt{qsx + Qix), that function is, in general, much 
narrower in the x direction than it is in the y direction. This fact causes 
the two-photon coincidence detection amplitude to not emulate the pump 
beam angular spectrum completely 47], |4^. For example, in the collinear 
{floe = fJ'eo = 0) degenerate monochromatic {ug = Ui = ujp/2) case, when the 
detectors are scanned in the same direction in the q space (q^ = = ^), 
the amplitude $00 is proportional to 



f (2q) sine {2ltq^) exp{-2iltq^) 



(61) 



It is clear that the angular spectrum will be clipped by the sine function and 
subjected to transverse walk-off. 

In type II phase matching, one of the down-converted fields has ordinary 
polarization and the other one has extraordinary polarization (e — )■ oe). In 
this case, one has two amplitudes: 
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It = apL/2, and l[ = (2ap — a)L/A ^ /t/2. The walk-off parameter a and 
the refractive indices n (ordinary) and fj (extraordinary) are calculated at the 
frequency uj = Up/2. Note the linear dependence of fioe and fieo on Us — uji- 
For the crossed cone configuration (ist . the state generated by type II SPDC 
is, according to Eq. (|63|) . 



where 



lV')sPDC ~ C'v \vac) + C2 |^)„ , (66) 

1^)11 ~ j j j j i^oe\(ls,^s,eo) \qi,UJi,ee) 

+$eo |q,,Ws,ee) |qi,Wi,eo)). (67) 

Again, if the crystal is thin enough and narrow-band filters are used, the 
state lip) II in collinear frequency-degenerate phase matching reduces to 



d<isdqi v{q, + qJdq,, e^) Iq^, eg) + |q„ e^) |q.. Go)). (68) 



D 



Eqs. fl62l) and (!63l) predict different transverse and longitudinal walk-off 
for ordinary and extraordinary down-converted photons. In type II phase 
matching, these effects can be considerable, and may decrease the quality of 
the entanglement between the photons. 

3.5. Coincidence detection probability and probability amplitude 

The coincidence count rate is proportional to the fourth-order correlation 
function 

C(p„ a) = (V^l E(-)(pJE(-)(p,)eW(pJeW(p,) I^) , (69) 

where E is the electric field operator and p^ and Pj are the transverse coor- 
dinates on the signal and idler detection planes respectively. Since lip) is a 



two-photon state, the correlation function fTMl) can be put in the form 41 

Cip„p,) = \^ip„p,)\^ (70) 

where C{pg,p^) can be interpreted as the two-photon detection probability 
and 

*(p„ p,) = (mc| eW(pJeW(p,) I^) , (71) 
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Figure 7: Illustration of an experiment to show the transfer of the angular spectrum from 
the pump beam to the spatial correlations of down converted photons J^] . 



as the two-photon detection amphtude. Eq. ( 1711) can be viewed as the two- 
photon, or bi-photon wavefunction. The electric field operator is given by 

E(+)(p) = C J dq, a(q) exp [^(q.p + ^/W^z)], (72) 

where C is a constant and z is the propagation distance. 

3. 6. The role of the spatial properties of the pump field 

The quantum states derived in section 13.41 show a dependence on the 
angular spectrum of the pump beam. This dependence indicates that the 
photon pairs can be prepared in a variety of quantum states through the 
manipulation of the pump beam. This preparation has measurable effects 
in the coincidence counting rate 16|, |49|. In order to illustrate this idea, we 



will describe one experiment in which the spatial shape of the coincidence 
distribution in the detection plane was prepared through manipulation of the 



pump beam |16| 



The set-up is shown in FigJTl The pump laser passes through a mask and 
a lens, before pumping the non-linear crystal. The mask has an aperture with 
the shape of the letter "C" . The lens is used to form the image of the mask 
in a plane situated after the crystal. The down-converted beams are detected 
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on the plane where the image of the mask is formed. In this condition, it is 
possible to recover the shape of the image "C" , in the transverse coincidence 
distribution. In order to image the coincidence distribution, the idler (or 
signal) detector is kept fixed, while the signal (or the idler) detector is scanned 
in two dimensions on the detection plane. A typical result is displayed in the 
inset in Fig. [71 The angular spectrum of the pump laser is transferred to the 
correlations between signal and idler photons, while the local intensities of 
each photon field are not affected 16l |. 

The coincidence count rate can be calculated using the state given in Eq. 
( 160|) together with several approximations simplifying the function $(q^,qj) 
given by Eq. (1571) . Assuming the thin crystal approximation, the SPDC 
state can be described by 



C2 dqjqi t;(q, + qj |q„ (X,)^ jq^, ai) 
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(73) 



where C2 is a constant. We neglect \vac) in Eq. ( !73|) . since it does not 
contribute to coincidence counts. 

Using Eqs. ([72]), ( 1^ and ( 172]) . the coincidence rate can be calculated: 
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(74) 



where >V(R; Zq) is the transverse spatial profile of the pump beam propa- 
gated to Z = Zq, and 



1 UJg 1 ^ UJi 1 

Zq Up Zs UJp Zi 

T, UJsZq Ui Zq 

R = Ps + Ps, 



(75) 



where Zg and Zi are distances between the crystal and signal and idler detec- 
tion planes, respectively. 

Eq. ( FM]) shows that the coincidence distribution depends on the Fres- 
nel propagator of the transverse profile of the pump field. This structure 
coincides with the image prepared in the pump, depending on the relation- 
ship between the wavelengths of the three fields and the free propagation 
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distances. This allows the engineering of spatial correlations by manipulat- 
ing the pump laser beam. Another important aspect is that the coincidence 
distribution depends on the sum of the signal and idler detector coordinates, 
and not either detector position alone, which is a signature of the spatial 
entanglement between the down-converted photons. 



4. Double slit experiments with twin photons 

Spatial properties of twin photons were experimentally investigated al- 
ready at the time of the first experimental observation of parametric down- 
conversion by Burnham and Weinberg [l|. They observed that the inten- 
sity correlations were stronger for certain combinations of detection angles. 
A more comprehensive study of these far field correlations were made by 
Grayson and Barbosa T^. Experiments exploring interference at a double 



slit experiment using twin photons and coincidence counting came in 1994 
[fill an d 1995 [3 and was further investigated in later work [s^ 34, 19, 17, 



52l . l53l . Il8l . l54l |55| . In this section, we review some of the key experiments 



exploring aspects of interference of photon pairs using a double slit. 



4.I. Nonlocal dependence of spatial coherence 

There are several possible ways to send one or two photons through a 
Young's double slit. Historically, the first choice is to send one of the photons- 
say the signal photon-through the double slit aperture, and let the idler 



propagate freely to the detection plane [l3|, as illustrated in Fig. [HI The 
two-photon coincidence distribution can be obtained in the same way as in 
Section 13.61 In order to take into account the presence of the slits in the 
signal beam, we include the diffraction integral in the electric field operator 
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EW(p) = J rfqc/q'a(q')r(q - q') exp 



^2 

q ■ p -Z -za 



(76) 



where T(q) is the Fourier transform of the double slit aperture A{p), za is 
the longitudinal coordinate of the slits, and Z = z — za- Using this operator 
and the quantum state (173!) in Eq. (I7T!) . we arrive at: 



dp'sdp'iMp's)^ ( Ip's + Ip'i, I e 
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Figure 8: Double slit experiment with twin photons. In Ref. [13|, the signal detector Dg 
was scanned to observe interference fringes. 



where W(p, za) is the transverse profile of the pump beam propagated to 
z = za and we have assumed that kg = ki = K/2 and Zg = Zi = z. Assuming 
that the distance Z is such that the Fraunhofer approximation is vahd, the 
detection probability is 



CiPs, Pi) oc 



k \ (k k 



(78) 



Approximating the aperture function by A{pg) = 6{xs + d) + S{xs — d) gives 
T(q) oc cos(g^(i). In this case the detection probability is 



C(P., Pi) oc 
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(79) 



From Eq. ( 1791) we see that the typical oscillations of a double slit interference 
pattern is found, as a function of the difference between the transverse coordi- 
nates of the signal and idler detectors. In Ref. [l3| , the coincidence patterns 
were obtained by scanning the signal detector Dg after the slits. A quali- 
tative comparison between the intensity fringes and the coincidence fringes 
observed as a function of the width of the idler detector aperture is shown 
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Figure 9: Qualitative comparison between intensity (single-photon) and coincidence fringes 
obtained by scanning the signal detector Ds as a function of the aperture of the idler 
detector Di. 



in Fig. O These results show that the visibihty of the coincidence fringes 
depends on the diameter of the idler detector, while the visibility of the in- 
tensity (single-photon) fringes depends only on the source properties. The 
visibility is high for a narrow aperture and low for a large aperture, implying 
a nonlocal dependence of the spatial coherence. The detection probability in 
this case is obtained from Eq. ( 1791) by integrating over the size of the idler 
detector 
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where D{pi) is the function describing the idler detector centered at po- 
sition cr in the detection plane. For a very narrow detection aperture, 
D{p^) ~ 5(pj — cr), and C£)(p^, pj) oc cos^[fc(i/22'(xs — a^;)], resulting in max- 
imum visibility. In other words, the post-selection by the narrow aperture in 
the idler mode projects the signal beam in a state with a narrower angular 
bandwidth. For a very broad detector, so that the width of D is much larger 
than the width of the angular spectrum v, we can approximate D ~ 1. In 
this case Eq. flHOj) gives a convolution, which washes out the oscillations of 
the cosine function, giving low visibility. A simple way of relating the idler 
detector diameter and the visibility of the coincidence fringes, following the 
lines of the classical van Cittert-Zernike theorem, was developed in Ref. 56 
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4-2. Ghost interference 

Another possibility in terms of double slit experiments using twin pho- 
tons and coincidence detection is the scanning of the idler detector instead 



of the signal detector IJ], as shown in Fig. [TOl The coincidence distribution 
displays interference fringes, while the intensity distribution of the idler field 
would never present interference fringes, since the double slit aperture is in 
the other beam. These fringes can be understood in terms of Eq. (17^ . which 
shows that the parameters of the coincidence interference pattern depend on 
the difference between the signal and idler detectors coordinates. Thus, by 
fixing the signal detector, and scanning the idler, the interference fringes ap- 
pear as a function of the position of the idler detector, even though the idler 
photon never passes through the double slit aperture. This type of exper- 
iment was named "ghost" interference, in reference to the "spooky" action 
at a distance attributed by Einstein to the non-local nature of entangled 
particles. 

4-2.1. The Advanced Wave Picture 

A method that aids in the understanding and design of experiments ex- 
ploring the spatial correlations between twin photons was introduced by A. 
V. Belinski and D. N. Klyshko \5l\. This method, known as the "Advanced 
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Figure 11: Advanced wave picture. 



Wave Picture" (AWP), consists of associating the coincidence detection and 
post selection to a temporal inversion of the propagation of one of the beams. 
In the case of the experiment of FigJTOl the AWP considers the system as if 
the fixed detector (signal) behind the slits was an incoherent light source and 
the scanning detector (idler) was registering the resulting intensity. The crys- 
tal plays the role of a mirror, as shown in Fig. [TT] The coincidence pattern 
obtained when the idler detector is scanned will be equal to the hypothetical 
intensity pattern obtained according to the AWP scheme described above. 
The AWP works independently of the optical components placed in the sig- 
nal and idler beams. It is important to notice that the crystal is treated as a 
mirror and the kind of mirror depends on the shape of the pump beam. This 



dependency was discussed in section 13.61 For example, in Ref. j49| it was 
shown that the curvature of the "mirror" is directly related to the curvature 
of the pump beam. 

4-3. The de Broglie wavelength of a two-photon wave packet 

Interesting possibilities arise when one sends both photons through the 
double slit aperture. Using the spatial correlations between the twin photons, 
it is possible to perform a double slit experiment in which the two photons 
always pass through the same slit 17|. In this way, the interference for a 
wave packet containing two photons can be observed. The signature of the 
two-photon interference is the spatial frequency of the fringes, which is twice 
the frequency of the fringes that arise in single photon interference. The 
spatial frequency is directly related to the de Broglie wavelength X^b = \/N 
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Figure 12: Experimental setup for two-photon interference at a double slit. The resulting 
interference pattern has oscillations corresponding to the spatial frequency of a two-photon 
wave packet. 



of the A^-photon wave packet [58]. An experiment which observed the de 



Broghe wavelength of a two-photon wave packet was reported in Ref. \V7 
and is sketched in Fig. [T2j A twin photon pair is produced in coUinear SPDC 
The beam containing the down-converted photons is sent through a double 
slit aperture with slit separation of 2d, after which the photons are split on 
a 50/50 beam splitter and detected in coincidence. In order to control the 
spatial correlations between the pair of photons in the plane of the slits, the 
pump beam is sent through a wire and a lens, placed before the crystal. The 
angular spectrum of the pump is produced so that the image of the wire is 
projected onto the double slit plane. The pump beam is blocked after the 
crystal, and does not actually pass through the slits. However the image of 
the wire is transferred to the correlations between the twin photons, such that 
there is a relative spatial separation 2d between the photons. Therefore, if 
one photon passes through one of the slits, the probability of finding its twin 
in the other slit is practically zero. Then, almost all photon pairs that pass 
through the slits must have passed through the same slit. The detection of the 
coincidence distribution is performed by scanning both detectors together, 
in order to mimic a two photon detector. Again, we can use the two-photon 
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Figure 13: Typical experimental results for two-photon interference at a double slit. The 
plots on the left show the profile of the pump beam at the double slit. The plots on the 
right show the two-photon interference pattern. In the lower pattern, the spatial frequency 
is doubled. 



state fl73|) and detection operator fl76|) for both the signal and idler fields to 



calculate the coincidence count distribution. We obtain 17 



C{x) oc \Bd{x)\^ + 4|5o(x)|2 + + (81) 

4Bd{x)Bo{x) cos + + 

\ Za Zl J 

45o(x)S_rf(x) cos + 

\za Zl J 



2Bd{x)B_d{x) cos 



^kx2d^^ 



where x is the position of the detectors, Bd, Bq and B^d are envelope func- 
tions which depend on the spatial amplitude distribution of the pump beam 
and the dimensions of the double slit aperture, za is the distance from the 
crystal to the double slit and zi is the distance from the slits to the detec- 
tors. In particular, Bj is proportional to the square root of the transverse 
profile of the pump beam at position x = j. We can see in Eq. ( 182|) that two 
oscillating terms depend on kx and one term depends upon 2kx, which is the 
two-photon interference term. Therefore, if Bq{x) = 0, then only the two- 
photon interference terms will survive. This condition is achieved through 
the preparation of the pump profile, using the wire and the lens to obtain a 
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zero intensity profile at x = 0. Then the coincidence rate will be given by: 



C{x) (X \Bd{x)\^ + \B_d{x)? + 2Bd{x)B^d{x) cos ( 1 • (§2) 

Fig. [T3] shows a computer simulation of typical experimental results, similar 



to those obtained in Ref. [17[. The type of two-photon interference fringes 



depends stro ngly upon the transverse profile of the pump beam at the double 



slit aperture [52[ . When the pump beam is a broad Gaussian distribution, the 
interference pattern depends only on the wavelength and not on the number 
of photons [59'], and an interference pattern, shown in the upper right-side 
of Fig. [131 is obtained. However, if one manipulates the pump beam so that 
the photons always pass through the same slit, the spatial frequency of the 
interference pattern changes to twice the original frequency, as illustrated in 
Fig. [T3b . Additional experiments reporting the observation of the de Broglie 



wavelength of two or more photons have been reported [60|, |61|, 162, |63|, |6J, |65 
The increase in spatial frequency in these experiments is equivalent to an 
increase in resolution, which may lead to interesting applications in optical 
lithography. 

4-3.1. Quantum Lithography 

We can conclude from the experiment discussed in section 14. 3[ that there 
are instances in which electromagnetic radiation with a given wavelength 
behaves as if it had a smaller wavelength, depending on the photon number. 
This idea has motivated the proposition of applications such as lithography- 
processes with increased resolution. It was called quantum lithography [66j, 



j m |68|, |69[ and is based on the possibility of creating two-photon wavepackets 
171]. In brief, using entangled pairs of photons and coincidence techniques, 
it is possible to obtain two-photon wavepackets which interfere and diffract 
in the same fashion as light with twice its wavelength. In principle, this idea 
can be extended to iV-photon wavepackets and in this case the increase in 
the resolution increases with A^. 

The term quantum lithography comes from the fact that the concept of 
photon is required to distinguish the behavior of one photon wavepackets and 
N photon wavepackets. However, there are propositions for the implementa- 



tion of quantum lithography using classical light sources [70|, |71|, 172|, |73!, [74 
In this case it is the interaction between the light and some material which 
would only absorb A^-photon wavepackets, and not J-photon wavepackets 
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(J 7^ A^). Detection of only A^-photon wavepackets projects the field onto 
the A^-photon component. 

Optical lithography is a process used to fabricate integrated circuits. It 
consists in the removal of a thin film from a substrate through the expo- 
sure to light. This results in drawings of very small circuits and electronic 
components. The future of integrated circuits in large scale depends on the 
possibility of reducing the size of these drawings. A limitation is given by 
the diffraction limit of the light used in the process, which depends upon the 
wavelength. Quantum lithography would allow one to obtain, for instance, 
the optical resolution of a UV beam, using near infra-red light. The main 
challenge towards quantum lithography seems to be finding proper materials 



which selectively absorb only A^-photon wavepackets [75 



An experiment demonstrating this idea, similar to the one described in 



section 14.31 was later reported in Ref. 69 



4.4- Non-local double slit 

Another possibility in terms of variations of a double slit experiment us- 



ing twin photons and coincidence detection is a non-local double slit [51 
illustrated in Fig. [TH In this experiment double-slit interference is observed, 
even though neither photon passes through a double-slit aperture. The non- 
local double slit is built from a single slit placed in the signal beam and a 
thin wire placed in the idler beam. The combined effect of these two com- 
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ponents results in double-slit type coincidence fringes. In order to calculate 
the coincidence count distribution, we can again use the two photon state 
given by Eq. f l73|) and the electric field operator given by Eq. f l76|) . where 
the transfer function T = 7i for the signal field corresponds to a single slit 
and the transfer function T = Ti for the idler corresponds to a thin wire. 
In order to observe the coincidence fringes for the non local double slit, it 
is necessary to focus the pump beam on the za plane. This allows one to 
approximate the pump beam amplitude transverse distribution by a delta 



function. The resulting coincidence count distribution is given by [51 



CiPs,Pi) oc 



J d^A.i^) A^i-^) exp (^tk^^ exp 



ZD - Za 



(83) 

where Ai and A2 are the transmission functions of the single slit and the 
thin wire respectively, za is the longitudinal coordinate of the single slit and 
thin wire, considered to be at equal distances from the crystal {z = 0), and 
zd is the position of the detection plane for both signal and idler detectors. 
The integral above is a Fresnel integral and the argument is given by the 
product of the functions Ai and A2 which results in the transmission function 
of a usual double slit aperture. Another characteristic which is present in 
this result, as well as many other SPDC experiments, is the dependence 
on the difference between the detector coordinates appearing in the Fresnel 
propagator in Eq. |83l 

4-5. The double-slit quantum eraser 

A double slit experiment with twin photons can be used to implement a 



quantum eraser [76|, |77[ , which shows the complementarity between wave-like 
and particle-like behavior. Exploiting the polarization correlation between 
down-converted photons, the idler photon can be used to herald either inter- 
ference fringes (wave-like) in the signal field or to determine which slit the 
signal photon passed through (particle-like). Furthermore, one can delay the 
decision about erasing this information or not. 

Typically, if one cannot determine through which slit a photon passes, 
then interference fringes will be observed in the intensity distribution after the 
double slit aperture. However, if one can mark the photon's path, so that the 
two paths are in principle distinguishable, then no interference fringes appear. 
It was previously argued that the disappearance of interference fringes was 
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Figure 15: The double slit quantum eraser for one photon. A/4 represents a quarter wave 
plate. 



due to the perturbation caused by the which-sht marker 78|. However, 



Scully, Englert and Walther showed that the interference fringes can reappear 



if one erases the which-slit information in the marker [77[. Hence the name 
quantum erasure. Let us begin by describing a simple double-slit quantum 
eraser with single photons. Let us suppose that a source emits single photons 
which propagate to a double slit. If interference fringes are visible on the 
detection screen after many photons have emitted, then one cannot determine 
through which slit each photon has passed. The state of the photons after 
the slits can be given by: 



71' 



tlphoton = -^im) + 1^2)], (84) 



where ipi and V'2 represent the passage through slit 1 and 2 respectively. In 
order to mark the path of each photon, let us place quarter wave-plates in 
front of each slit, so that if the horizontally polarized photon crosses one 
of the slits, it leaves circularly polarized to the right and if it crosses the 
other it leaves circularly polarized to the left, as in FigJTSk. Then, a mea- 
surement of the photon polarization is enough to determine through which 
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slit it has passed. This causes the interference pattern to disappear, even if 
the polarization is not measured. Even the possibility of obtaining which- 
path information is enough to destroy the interference, since the presence of 
the wave plates in each slit entangles states ipi and 1^2 with the polarization 
degree of freedom, so that the state becomes: 



1 

3^ 



)lphoton = -^[\\R)\iIJi) + \L)\^2)], (85) 



where \R) and \L) are right and left circular polarization states respectively. 
The which-path information is available due to the entanglement between 
the state describing the passage through each slit and the output polariza- 
tion state. It can be erased by projecting the two orthogonal polarization 
states l-R) and \L) onto a-say-linear polarization state, | if) (horizontal) or 
I V^) (vertical) for instance, as in Fig. [T5b . Therefore, after the polarizer, the 
which-path information is no longer available, and the interference fringes 
are recovered. 

Using pairs of photons allows for a delayed choice quantum eraser, as 



reported in Refs. |53l . ISJ, l79||. A sketch of an experiment is shown in Fig. [161 
A pair of photons is prepared in a polarization-entangled state: 

1 
71 



Q^ = —[\H),\V)s + \V),\H)4. (86) 



The signal photon is sent through the double slit aperture and quarter-wave 
plates. The state of the photons after passage through the slits is given by: 

\ij)QE = ^mmslA) + \R)s\i^2)) + t\VU\R)s\A) + \L)s\i^2))]. (87) 

The polarization entanglement allows the idler photon to be used as a which- 
path marker or eraser. If it is projected onto polarization H or V the signal 
photon will be projected onto state \L)s\ipi) + 1-^)81^^2) or \R)s\4'i) + \L)s\4'2) 
respectively. Therefore, no interference is observed. However, if one projects 
the idler photon onto |-|-)(+45 degrees) or |— )(-45 degrees) linear polariza- 
tions, the which-path information is erased and the interference fringes re- 
turn, as illustrated in Fig. [T71 In order to see this, we write the state ( 1871) 
in the +/— basis: 



1 



)qe = — [|+),,|+),(|^i) - ^1^2)) + ^l-).l-)s(l^i) + (88) 
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Figure 16: The delayed choice double slit quantum eraser. Quarter wave-plates are used 
to mark the path of the signal photon. Which-path information can be recovered or erased 
by projecting the idler photon onto the appropriate polarization state. 



We see in Eq. ( I88l) that projections of the idler photon onto or |— ), 
projects the signal photon onto a superposition of {ipi) and |'02) states, which 
gives interference fringes. This quantum eraser allows for "delayed choice", 
meaning that the decision about the measurement of the idler photon can be 
delayed until the signal has already been detected. 

5. Quantum imaging 

In sections 13.61 and we discussed experiments in which the image or in- 
terference pattern of an object appears in the two-photon coincidence counts. 
These types of experiments have led to the idea of "quantum imaging" , which 
in some cases may present advantages over classical imaging techniques, such 



as increased spatial resolution [80|, ImI, |82|, |83| , as discussed in section 1^?^ The 



term quantum images generally refers to images that appear in two-photon 
coincidence distributions. They are called quantum, because they were first 
produced with the quantum-correlated twin photons produced in parametric 
down-conversion. However, it is also possible to produce correlated images 
with classically correlated light sources. In this section we will discuss quan- 
tum and classically correlated images and some proposed applications. 
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Figure 17: The delayed choice double slit quantum eraser. Each polarization projection 
of the idler photon corresponds to a state preparation of the signal photon. Projections 
onto vertical and horizontal polarizations results in which-path information. Projections 
onto diagonal polarizations results in erasure of which-path information and interference. 



5.1. Quantum images 

The image of a double-slit aperture can be understood in terms of the 
diffraction theory of classical optics as the near- field distribution of the double 
slit, whereas the interference pattern is related to the far- field distribution. 
The amplitude distributions of the electric field in these two cases are related 
by a Fourier transform. This is also true for the correlated images and in- 
terference patterns which may appear in coincidence distributions. However, 
the propagation and detection of both signal and idler fields may influence 
the resulting spatial distribution. 

In order to introduce the idea of a quantum image, let us consider the 
experiment reported by Pittman et al. 15| and illustrated in FigJTSl Para- 
metric down-conversion is produced as usual. An aperture mask is placed 
in the signal beam close to the detector and a lens also placed in the signal 
beam. Coincidence detection is performed such that the signal detector is 
a fixed large-aperture bucket detector, while the idler detector is scanned in 
the transverse plane. A 3-D plot of the coincidence counts as a function 
of the transverse coordinates x and y reproduces the image of the aperture 
ABC. See the inset in the upper part of Fig. [181 The lens in the signal 
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Figure 18: Observation of a quantum image, similar to the experiment of Pittman et al. 
[iSj . The image of the ABC mask appears in the coincidence counts when the detector 
Dl is scanned. 



beam is placed according to the advanced wave picture (see subsection 14.21 
for details). The imaging condition given by the thin lens formula j = 7 + ~ 
is applied, where / is the focal length, o is the distance between the lens and 
the aperture mask (the image plane) and i is the distance between the lens 
and the crystal plus the distance between the crystal and the idler detector. 
According to the advanced wave picture, the crystal plays the role of a mirror 
in the imaging process. 

There are some particular aspects which distinguish this image pattern 
from ordinary intensity images. One feature is that the image is observed in 
coincidence as a function of the position of the idler detector, even though 
the idler did not pass through the aperture mask. The image does not appear 
in the intensity distribution of the idler alone, but only when it is counted in 
coincidence with a detection of the signal photon. The intensity distribution 
depends only on the source shape as usual, as illustrated in the inset in 
the lower part of Fig. [181 Another aspect is that the dimensions and the 
resolution of the quantum image depend on the pump beam wavelength and 
spatial properties. 

Initially, it was widely believed that the quantum correlation of the down- 
converted photons was necessary to observe quantum images 8J|. Thus 



it 



was quite surprising when similar types of correlated images were produced 
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Figure 19: Observation of a classically correlated image using a laser pulse as reported in 
Ref. 



with classical sources. This is the subject of the next section. 

5.2. Classical two-photon images 

A g reat deal of discussion was triggered by the work of Bennink et al. 



85l . l86[ , in which it was demonstrated that it is possible to obtain correlated 
images using classical light, in the same fashion as it is done with entangled 
photons. Bennink et al.'s experiment is sketched in FigJTOl A laser pulse is 
sent to a mirror and then to a 50-50 beam splitter. The transmitted beam 
propagates through an aperture mask and is detected. Let us call this the 
signal beam. The idler beam is reflected at the beam splitter and sent directly 
to a CCD camera, which records the spatial intensity distribution. It is not 
necessary to use a pulsed laser, the pulse-like scheme can be obtained for 
instance with a chopper in a c.w. laser. However, it is important to have 
pulses, so that there will be a temporal correlation between signal and idler 
pulses. 

The angle of the mirror before the beam splitter is randomly modulated 
in both axial and azimuthal directions. This artificially broadens the angular 
spectrum of the light sent to the beam splitter, and the temporal correlation 
results in an angular correlation between the signal and idler beams. The sig- 
nal detector after the aperture mask is kept fixed, while the idler is detected 
with the CCD camera which only records when a trigger pulse is sent by the 
signal. The triggered CCD image reproduces the ABC aperture, just as in 
the case of correlated photon pairs from SPDC. Almost all ingredients natu- 
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Figure 20: Observation of a classically correlated image using a thermal light source. 



rally present in the pair of photons produced in parametric down- conversion, 
are artificially introduced here: i) temporal correlation and ii) angular corre- 
lation. However, since the correlated pulses are produced classically, there is 
no entanglement. This experiment demonstrated that it is possible to obtain 
a correlated image using classical light, and opened up the discussion about 
the role of entanglement in quantum imaging experiments. 

The production of correlated irnages with classical light was extended to 



thermal/chaotic light sources 0, H H H 
correlated images were not actually quantum 



90 



supporting the idea that 
A typical experimental set- 
up for the observation of correlated images with thermal light is shown In 
Fig. [20] and is similar to Fig. [191 However, the pulsing to obtain temporal 
correlation and spatial modulation to obtain spatial correlation is not neces- 
sary anymore. Due to the natural bunching of photons in thermal light, it is 
possible to obtain temporal and spatial correlation directly from the thermal 
light. This kind of correlation was already observed in the famous exper- 
iment by Hanbury Brown and Twiss (HBT)[91]. In the HBT experiment, 
sketched in Fig. [HI the light is sent to a beam splitter and intensity corre- 
lations are measured. One detector is kept fixed, while the other is scanned 
along the detection plane. Results show that the width of the coincidence 
distribution is proportional to the spatial coherence of the input field. The 
larger the coincidence distribution, the smaller the spatial correlation and 
the larger the spatial coherence. For incoherent sources like thermal sources, 
the spatial coherence depends only on the source dimensions, as discussed in 
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Figure 21: Sketch of the Hanbury Brown and Twiss experiment. C is the intensity corre- 
lation. It is proportional to fii2, the normalized spatial degree of coherence. 




section 13.1.21 In this way the HBT interferometer was used to measure the 
diameter of stars. Correlated imaging with a laser source required a modu- 
lation to induce spatial correlations since the spatial correlation is inversely 
proportional to the spatial coherence, which in the case of a laser source can 
be quite large. 

While classical light can also be used to produce correlated images, there 
is a price to be paid: the signal to noise ratio is generally higher than for 



the output state of a down-converter [87|, [92] • Quantum imaging has been 



analyzed in the general case of light fields with Gaussian photon statistics 



93l . |94| , and this formalism has been applied to two-photon imaging [95 



5.3. Phase objects 

The discussion about the differences between quantum and classically 
correlated images has led to an important test, proposed and realized exper- 
imentally by Abouraddy et al.[96]. The authors observed correlated images 
using phase objects. Fig. [22] shows a sketch of a typical experimental set-up. 
The transmission configuration using an aperture is replaced by a reflection 
configuration using a micro electromechanical system (MEMS) micro-mirror. 
The spatial amplitude distribution of the reflected beam is not changed, while 
the spatial phase distribution is modulated. The signal photon propagates 
through this plate and is detected afterwards by a flxed detector Dl, outfltted 
with a small pinhole. The idler photon propagates directly to detector D2, 
which is scanned in the transverse plane. The spatial coincidence distribution 
will be determined by the diffraction pattern of the signal beam. A typical 
coincidence proflle is illustrated in the inset of Fig. [25J This kind of two 
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Figure 22: Observation of a quantum image with a pure phase object. 



peak profile, is obtained when a double-slit object with phase vr separated 
by a line of phase is prepared in the MEMS. It was shown in Ref.|93], that 
coherent imaging can also be achieved with pseudo-thermal sources. The 
observation of quantum images of phase objects is a further demonstration 
of the differences between the correlated images with classical and quantum 
light. 

5.4- Spatial resolution of magnified images: quantum versus classical 

It is well know that the spatial resolution of an image increases when the 
wavelength of light used for illuminating the object is decreased. On the 
other hand, short wavelengths in the UV range or smaller usually damage 
an illuminated object. Needless to say, it would be advantageous to achieve 
the spatial resolution of short wavelengths while illuminating the object with 
a light beam with longer wavelength. Spatially entangled two-photon light 
beams generated by SPDC are good candidates for this special light source 
because it has two wavelengths associated with the two-photon beams: the 
central wavelengths of the individual photons A and the de Broglie wavelength 
associated with the biphoton wavepacket equal to A/2 , as discussed in section 
14. 3[ It was demonstrated experimentally in Ref. [Slj that magnified images 



of objects illuminated by two-photon wavepackets, as in a microscope [98 



can be generated with spatial resolution better than the diffraction limit and 



are self-apodized [99 



For a comparison with a classical imaging experiment suppose an object 
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is illuminated by a classical coherent light source. The image is produced 
with a lens with a focal length /, separated from the object by a distance 
zl. The image plane is located at a distance zn far from the lens. If it is 
assumed that the light arriving at the object is a plane wave, the electric 
field at the image plane is 



\Zl Zd) 



E{p) = OmJ^ + ^)dt (89) 



where O(^) is the transmission function of the object, p is the transverse 
position in the image plane, k = 27r/A is the wave vector of the incident 
plane wave and Tl is the Fourier transform of the magnitude of the lens 
transmission function. This last term limits the spatial resolution of the 



image due to the finite size of the lens aperture [99 . 

In the imaging experiment with photon pairs generated by SPDC, the 
object information appears in the two-photon probability amplitude at the 



image plane. It was shown theoretically in Ref. [80| that the probability 
amplitude of simultaneously detecting two down-converted photons in the 
image plane of an object illuminated by the signal photon is 

np) = / ^(^)^^ + ^) (90) 

where it is assumed that zl = zu = zls, zd = z^i = z^g. Vector p describes 
the detection position at both signal and idler detection planes, so that Pj = 
—Pg = p. This last assumption means that image is measured by displacing 
the signal and idler detectors simultaneously in opposite directions. It has 
also been assumed that the pump laser beam at the object plane has a 
transverse profile which is sufficiently narrow to be described by a Dirac 
delta function. Tp is the Fourier transform of the function F defined as 

F(v) = j |Aii(u + v)Pi,(u-v)|du, (91) 

where \Als\ and {Anl are the magnitude of the aperture functions of the 
lenses placed in signal and idler paths, respectively. F{\) is the correlation 
of the lens transmission function magnitudes [soj. In comparing Eqs. f l89|) 
and fl90|) . one notices the presence of the factor 2k in Eq. fl90|) instead of k as 
in Eq. f l89|) . Eq. fl90|l is equivalent to an image generated by single photons, 



however with wavelength equal to the De Broglie wavelength [58|, |17[ of the 
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biphoton A/2, and thus shows improvement in the spatial resolution. On the 
other hand, the function F present in Eq. (M?]) describes an effective lens 
with aperture transmission function that has a magnitude that is equal to 
the correlation of the magnitudes of the aperture transmission functions of 
the actual lenses, as it can be seen in Eq. (|9T]) . The function F describes an 



effective apodized lens (99 



The image generated by the parametric down-converted photons is better 
resolved than a similar one generated by an infrared classical light source with 
the same wavelength, and is resolved better than the diffraction limit in this 
case. In spite of this, the quantum image resolution is not as good as the 
resolution of the image produced by the pump laser source (Ap = A/2), since 
the transmission function of the effective lens seen by the photon pairs is not 
equal to the transmission function of the original lenses. Apodization effects 
are also observed in the twin photon image without physically apodizing the 
lenses used, which could lead to interesting applications. 

The spatial resolution of images obtained in quantum fourth-order imag- 
ing has also been compared with that obtained in a classical second-order 



incoherent imaging method [82[. The intensity at the image plane of an 



object illuminated by a classical incoherent light source is [36 



/(P)= + (92) 

where Tp is the Fourier transform of the function p, defined as 

p(u) = y"AL(u + v)AL(u- v)c/v, (93) 

= 27r/A, in which A is the central wavelength of the incoherent light beam 
and Al'is the aperture transmission function of the lens. The function p(u) 
is the auto-correlation function of the aperture function of the lens used for 
imaging, and describes the aperture function of an effective lens. Note that 
in the two-photon image the effective lens transmission function f l9ip is not 
an auto-correlation function but a cross correlation function of the signal 
and idler aperture lens transmission functions. Another difference with re- 
spect to the object function is that in the incoherent imaging only the square 
modulus of the object function appears in the optical intensity expression, 
i.e., this method is not sensible to the object phase. In quantum imaging 
the complete object function is present in the coincidence distribution at the 
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image plane, i.e., the image carries the information about the object spatial 
phase. A situation where the object function has no phase variation was an- 



alyzed in Ref. [82|, but the two imaging methods still produced completely 
different results. The transmission function of the lens pair in the two-photon 
imaging was modified such that the cross- correlation function between them 
was different from the auto-correlation of one of the lenses. When the object 
was illuminated by an incoherent classical light source, the image was not 
resolved spatially. On the other hand, an image with very high spatial reso- 
lution was measured when the object was illuminated by one down-converted 
photon, and both photons were detected in coincidence at the image plane. 

6. Spatial correlations: quantum versus classical 

In the 1990 's, many novel and interesting features of spatial correlations 
were observed experimentally and well described by the quantum theory. 
Many of these results suggested that the spatial correlations were indeed 
quantum correlations, even though the non-existence of a classical analog 
was not discussed in most of these contributions. For example, the non-local 



conditionality of two-photon interference fringes |14j . |60| was well described 
by quantum mechanics. However, there was never any formal proof that 
the conditional interference fringes could not be exactly reproduced by some 
special classical light source. In Ref. jsj, it was shown that the simple detec- 
tion of coincidence counts above the accidental coincidence rate is a proof of 
non-classical behavior. However, this refers to the temporal correlations and 
cannot be directly extended to spatial variables. 

Perhaps the first experimental demonstration of the quantum nature of 
light via a fourth-order correlation function was photon ant i-bunching, first 



observed in the time domain by Kimble et al. in 1977 |100| . Demonstration of 
the non-classical nature of photon anti-bunching was based on the violation 
of a classical inequality. In the same spirit, the anti-bunching of photons 



was later demonstrated in the transverse spatial domain in Ref. [18|, and 
was probably the first formal demonstration of the non-classical behavior 
of spatial correlations. In 2004, Howell et al. 2lj| and D'Angelo et al. 22 



demonstrated that the pairs of down-converted photons are indeed entangled 
by experimentally violating classical separability criteria. 

In this section we review several of these important results concerning the 
non-classical nature of transverse correlations. The observation, detection 
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Figure 23: Double-slit interference with two photons. 



and quantification of spatial entanglement will be discussed in sections [7] and 

El 

6.1. Conditional interference and complementarity in one and two-particle 
interference patterns 

An interesting idealized two-particle double-slit interference experiment 



was discussed by Greenberger, Horne and Zeilinger |101l | in 1993. In this 
gedanken experiment, two identical double-slits were placed in opposite sides 
of a linear source Vt with extension I that always emits two particles with mo- 
menta approximately equal and opposite, as sketched in Fig. |23l The angle 
between the slits and the source is 6 and, due to the momentum correlation, 
the particles pass either through the slits A and A' or B and B' . The arrival 
of the particles is registered at two opposite screens after each double-slit. 
For the case where the daughter particles are photons and / >> A/6', the 
probability am plitu de for one of the particles to arrive at P and the other to 
arrive at P' is |l01 



\l'2(xj, Xs) oc cos[k6\xi — Xs)] (94) 

where Xs and Xi are the transverse coordinates, /c = 27r/A, A is the wavelength 
(de Broghe or optical) related to the emitted particles, and 6' is the angle 



50 




Figure 24: Double-slit interference with two photons and parametric down-conversion. 
Focusing the pump beam with a lens induces an anti-correlation (dependence on the sum 
of coordinates) and using a a lens plus a wire induces correlation (dependence on the 
difference of coordinates) . 



that is subtended by the hole pairs and the detecting screens. On the other 
hand, if / << A/6', the probabihty amphtude is 

'^2{xi,Xs) oc cos{k9' Xi) X cos{k9' Xs). (95) 

Expression shows that the coincidence interference pattern detected 
when / >> X/6 presents "conditional fringes". If the interference pattern 
is recorded by two movable detectors Z)j and Ds, the term "Conditionality" 
means that the position of the interference fringes will vary depending on 
the position of both detectors. Suppose the interference pattern is recorded 
by keeping one of the detectors (-Dj, for example) fixed at some position x 
while Ds is scanned. For two different positions x, two displaced interference 
pattern will be recorded. Expression (l95il shows that for I « X/9 a. product 
of independent single-particle interference patterns ("independent fringes") 
is measured and conditional fringes are not detected. This condition is the 
requirement to detect the usual Young single-particle interference pattern. 

us- 



This idealized experiment was tested experimentally in reference [52 



ing photons obtained from SPDC The experiment is illustrated in Fig. 
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Two-photon interfer ence experim ents with spatial interference patterns were 
also observed in 0, Il02 . 14, 56, 103 . 17 1 but as discussed above 101 1 the 
detection of a fourth-order interference pattern does not guarantee the pres- 
ence of conditional fringes. Using the quantum multimode theory presented 
in section [3] to calculate the fourth-order correlation function as a function 
of the detector positions, the number of coincident photons at positions Xi 



and is 52 



Nc(xi,Xs) oc A{xi,Xs) + 2 Bi{xi,Xs) B2{xi,Xs) cos 

+2Bi{Xi,Xs) B3{Xi,Xs) COS 
+2B2{Xi,Xs) B^{Xi,Xs) COS 

+2B3{xi,Xs) Bi{xi,Xs) cos 

+2B2{Xi,Xs) B3{Xi,Xs) COS 

+2 Bi(xi,Xs) BJxi,Xs) cos 
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with 



A{xi 



Xq 



B\ {Xi^ X g 



B2{Xi,Xs 



B3{xi,Xs) \^ + \B^{xi,Xs) 1^(97) 



Here 2d is the separation of the double-slits, 2a is the width of each slit. The 
first four cross terms in f l97|) give independent fringes that depend only on Xs 
or Xi, while the last two terms give conditional fringes which depend upon 
± Xj. The interesting fact that appears from this calculation is that the 
diffraction terms Bj{xi,Xs), j = 1,2,3,4, are proportional to the transverse 
profile of the pump laser at different transverse positions in the plane of 



2,3, Bi{xi, 



X a 



OC 



the double slit aperture: Bj{xi,Xs) oc }V{0,za) with j 
yV{d,ZA) and B^^XijXg) oc yV{—d,ZA)- Generating the photon pairs using 
pump beams with different transverse profiles at the double-slit plane, the 
authors in Ref. ^2| were able to measure two-particle independent fringes, 
as well as conditional fringes that depend on the difference or the sum of the 
detector position coordinates. Conditional fringes which depend on the sum 



52 



of the po sition coordinates have been predicted theoretically in references 



104, 105 



Greenberger, Horne and Zeilinger 10 11 ]. analyzing the two extreme cases 
discussed above (independent versus conditional fringes) affirmed: "... there 
is a sort of complementarity between one- and two-particle fringes: the con- 
dition for seeing one precludes the possibility of seeing the other". When 
the source dimension or the distance from the double-slits to the source is 
such that the idealized experiment is analyzed in an intermediate geometry 
when compared with the extreme cases mentioned above, both single-particl e 



fringes and conditional two-particle fringes are present |106l . llOll Il05l |107 



Horne |105| derived a complementarity relation for the one- (vi) and the 
two- particle (f 12) fringe visibility for a general state of a two-particle system 
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V, + V 



12 



(98) 



The single-particle fringe visibility is derived from the normalized single- 
particle probability density, which is obtained by integrating the joint prob- 
ability density with respect to the position of one of the particles. The 
two-particle visibility is obtained from the two-particle "corrected" probabil- 
ity density 1091 . Il08l |. This definition prevents the visibility Vu from being 
unity when partic les 1 and 2 are prepared in a product state. Abouraddy 
and collaborators llOl . Illl| studied the two-particle double-slit interferome- 



ter for the case where the particles are generated by means of SPDC. Photon 
pairs generated collinearly by the crystal are passed through a double-slit 
aperture and are imaged by a 2/ lens system at the coincidence detector 
camera. For this geome try, t he authors verified the above complementarity 

Calculati on o f vi and V12 from the measured 
The single-photon coincidence 



relation experimentally |110 



data follows the strategy used by Horne |105 



probability associated with the two-photon pair transmitted by the double- 
slit is calculated by integrating the fourth-order correlation functions with 
respect to one of the position variables. 

6.2. Spatial antibunching 

The experimental observation of the spatial anti-bunching of photons con- 
sists in the production of a homogeneous optical field for which it is more 
probable to detect coincident photons at spatially-separated positi ons i n the 
transverse plane. The original idea of anti-bunching of photons 100l | was 
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connected to the time interval between the emission of two photons by a 
hght source. Spatial anti-bunching means that the probability of finding two 
photons together at the same position in the transverse propagation plane 
is smaller than the probability of finding them at different positions. For a 
classical light beam, the best that one can do is to have equal probability of 
detecting photons together or apart. A thermal source produces light that 
is bunched in time and in space, while an ideal laser produces light which 
is neither bunched nor anti-bunched. It is also important to note that both 
bunching and anti-bunching require that the fields be stationary (depend 
only on r = ti — and homogenous (depend only on 6 = Pi — ^2)- If this 
were not the case, it would be possible to simulate temporal anti-bunching 
using a low intensity pulsed laser, for instance. 

From a mathematical point of view, the inequality 

r(2'2)(5,r) < r(2'2)(0,0), (99) 

must be satisfied by any classical field. F*^^'^-* is the fourth-order coherence 
function, which is proportional to the coincidence count probability. The 
violation of this inequality implies anti-bunching if the condition of station- 
arity and homogeneity is fulfilled. This means that the temporal and spatial 
arguments of all coherence functions {t,t + r, p, p + S) can be replaced by 



(r, 5). The experiment performed in Ref. [18| is sketched in Fig. The 
non-linear crystal is pumped by a Gaussian profile laser beam and collinear 
type II down-conversion is produced. Signal and idler photons are sent to 
a double-slit aperture equipped with quarter wave-plates (see below). After 
propagation through the slits, the photons are sent to a 50/50 beam split- 
ter (BS) and single-photon detectors are placed in each output port. The 
detectors are carefully calibrated, so that the transverse coordinates in both 
detection planes are equivalent. This arrangement is essentially a two-photon 
detector. Detector Dl is kept fixed while D2 is scanned in the transverse 
plane. For an ordinary double-slit, the resulting coincidence pattern would 
display the usual interference fringes, equivalent to the intensity fringes in a 
double-slit experiment. To produce a spatially anti-bunched field, a phase 
shift of vr is produced in the coincidence fringes, so that the interference max- 
ima are now minima and vice versa. In this case, the probability of finding 
two photons together in the same point {6 = 0) is smaller than the proba- 
bility of finding them apart. This is the signature of spatial anti-bunching of 
photons, provided the condition of homogeneity is fulfilled. 
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Figure 25; Experiment for the observation of the spatial anti-bunching. The vertical axes 
in the plots are normalized coincidence rates for the solid lines and single count rates for 
the dashed lines. The horizontal axes are D2 detector positions in arbitrary units. The 
position of detector Dl is fixed in A) -2.5, B) and C) -1-2.5 and D2 is scanned. 



The 71 phase shift is obtained using a double-sht aperture with a zero- 
order quarter- wave plate placed in front of each slit. The fast axis of one 
wave plate is oriented along the vertical direction and the fast axis of the 
other is oriented along the horizontal direction. The pump beam is focused 
in the plane of the slits, which causes an anti-correlation in the positions of 
the photons in this plane 16|, so that if one photon crosses one slit, its twin 
will pass through the other slit. The type-II SPDC produces photons in the 
polarization state \H) ^\V). The position anti-correlation at the plane of the 
double-slit guarantees that the photons pass through different slits, resulting 
in the state 



1^) 



1 



AB 



-^(IV'i, H)\i^,, V) + e^^/V2, H)e'^/'\ij,, V)), 



(100) 



where ipi and are the quantum states describing the spatial variables when 
the photons pass through slit 1 or slit 2, respectively. The it/ 2 phase factors 
are due to the delay introduced by the zero-order quarter-wave plates. The 
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coincidence probability is [l8|, [19 

C(Xs, Xj) OC 1 — COS 



kd 



{Xi - X2) 



(101) 



where k is the wave number of the down-converted fields, d is the separa- 
tion between the slits, z is the propagation distance, and xi and X2 are the 
positions of detectors Dl and D2. The coincidence probability is zero when 
Xi = X2, corresponding to a spatially anti-bunched field. Typical coincidence 
and intensity measurements observed in Ref. 18| are shown in the insets of 
Fig. [25l The coincidence distribution is near zero at when both detectors 
observe the same point, and the single-photon intensity is nearly constant 
for all positions. The coincidence distributions violate inequality (!99|l . as 
r'^^'^^(<5) > r(^'^)(0), for many values of 6. The measurements are performed 
for several positions in the transverse planes of both detectors, demonstrating 
homogeneity. It was also demonstrated that the same effect can be observed 
without the use of a double-slit, with free propagating signal and idler beams 
Spatial anti-bunching has also been produced using Hong- Ou-M andel 

This 
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interference of photons in an anti-symmetric polarization state 
arrangement allows for the production of a beam of photons in a 
polarization state. 
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■'singlet" 



7. Spatial Entanglement 



Quantum entanglement is a correlation that can exist between two or 
more quantum systems. In the case of pure states, entanglement implies that 
the quantum state 1^1') ^2 describing the-say-two systems cannot be separated 
into a product of quantum states lip) ^ \(f)) 2- If the bipartite state can be 
written in this form it is said to be separable. In the more general case of 
mixed states, a separable bipartite state can be wr itten as a convex sum of 



tensor products of local density matrices |114j . |115 



^sep 



(102) 



where YliPi ~ ^ Pi — 0- ^ bipartite quantum state cannot be written 
in the form f ll02p . it is entangled. For a more in depth discussion of quantum 
entanglem ent in general, we refer the reader to two recent review articles 



114, 115 



56 



As discussed in sections IH El although many experiments performed 
in the 1990 's strongly suggested that photon pairs produced by SPDC were 
indeed entangled in the spatial degrees of freedom, it wasn't until 2004 that 
experimental tests confirmed this fact 2l|, |22| . These experiments made use 



of entanglement criteria involving meas urements of the sum and difference of 
the position and momentu m variab les 
conditional measurements 



118, 119 



1161 . [llTj, or EPR criteria involving 



7.1. Continuous Variable Entanglement Criteria 

There has been much work concerning the detection and characterization 
of entanglement in continuous variable degrees of freedom. An overview of 
the continuous variable formalism, and entanglement criteria f or continuous 



variables, can be found in a number of review articles |l20l . Il2lj . Though the 
physical system of interest is usually the field quadratures of intense fields, 
continuous variable entanglement conditions are also applicable to spatial 
entanglement of photon pairs. 

For the spatial degrees of freedom of photon pairs, it is most useful to 
consider the global operators 

X-t = Xi ± X2, (103a) 

P± = Pi ± P2, (103b) 

where in the present case x and p are the position and momentum operators 
for photons 1 and 2, and we assume that [xj, p;.] = idj^k and j, k = 1,2. In the 
case of spatially entangled photons, the position and momentum observables 
correspond to measurements in the near and in the far-field respectively, 
relative to the source plane, namely the SPDC crystal. The near field is 
associated with a position measurement at the source plane, while the far 
field is associated with a momentum measurement at the source plane. 

In most cases, entanglement can be identified through violation of one 
of a number of inequalities, which thus serve as entanglement witnesses. 
For instance, the Mancini-Giovannetti-Vitali-Tombesi (MGVT) criteria [117| 
reads 

(A2x±)(A2p^) > 1, (104) 
where stands for the variance, while the separability criteria of Duan, 



Giedke, Cirac and Zoller (DGCZ) [116 



IS 



(105) 
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where a is a local scaling parameter which guarantees that the quantities 
in the sum are dimensionless. One can opti mize condition fllUSp over a, 
in which case one arrives at condition f ll04p [l22'|. All separable states of 
the form fll02p will obey inequalities fll04p and fllOSp . and thus violation of 
either of them is a sufficient condition for identifying quantum entanglement. 
Inequalities fll04p . f llOSp are examples of a rnore general class of entanglement 
witnesses involving second-order moments 123| , for which Hyllus and Eisert 
have provided an optimi zatio n procedure 122| |. 

The Simon criterion jl24| is a necessary and sufficient condition for Ix N 
mode Gaussian states [125|. To evaluate this criterion, it is necessary to 



reconstruct the covariance matrix 120, 121 . Violation of the Simon, MGVT 



or DGCZ criteria are suf ficien t for i dent ifying quantum states with negative 
partial transpose (NPT)[i26|, Il27l . |128|. In the general case, violation of 



these criteria is not a necessary condition for entanglement even in the case 
of Gaussian states, due to the existence of bound entangled states which have 
a positive partial transpose. 

For non-Gaussian states, these second-order criteria are sufficient, but 
not necessary, and there exist NPT entangled states which do not violate any 
second-order criteria. In this case one can use the NPT criteria of Shchukin 
and Vogel, which provides a hierarch y of i nequalities involving combinations 
of second and higher-order moments 129l |. In fact, they have shown that in- 
equalities f ll04p . fllOSp and the Simon criterion are special cases of this general 
NPT criteria. Violation of any inequality in the Shchukin- Vogel hierarchy 
indicates that the quantum state is NPT, and is thus a sufficient condition for 
entanglement. Higher-order inequalities obtained from the Shchukin- Vogel 
criterion have been violated recently for a spatially non- Gaus sian two-photon 
state which does not violate any second-order criteria 130|. We note that 



there exis t several other CV entanglemen t criteria invo l ving higher order 



moments 



i3iL I132L liaaL UM uM uM uM UM uM uM Ml 



7.2. Correlations in the near and far field 

The CV entanglement criteria discussed in section 17.11 can be easily ap- 
plied to the spatial correlations of SPDC photon pairs. In this case, the posi- 
tion variable x corresponds to the transverse position p and the momentum 
variable p is related to the transverse wave vector via p = hq. The position 
and wave vector distributions can be observed through intensity measure- 
ments in the near and far field of the non-linear crystal, respectively. Using 
simple optical systems composed of lenses and free space, it is possible to 
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measure the spatial distributions of the down-converted photons in the near 
and far field. In the near field, one generally observes a strong correlation 
between the positions of the signal and idler photons. An intuitive way to 
understand this correlation is through the fact that the two down-converted 
photons are "born" from the same pump photon simultaneously. Thus, the 
photons are produced at approximately the same transverse position inside 
the crystal. If the crystal were infinitesimally thin in the longitudinal direc- 
tion, this would produce a point-like correlation between their positions in 
the crystal plane. 

The field operator which describes detection in the near-field is given by 

EM oc J dpaiq)e-"'f. (106) 



Assuming that the quantum state of SPDC photon pairs is given by Eq. (!59|) . 
the near-field correlations are governed by the profile of the pump beam 
>V(p) and the Fourier transform of the phase matching function: T{p) = 
•^{i{q)}- Typically, r(p) is much narrower than W(p) at the crystal face. 
The detection amplitude is defined as \E'nf(p^, = {0\E^ {p^)E^ (p^)]'^) and 
its square modulus is proportional to the coincidence rate. 

If W(p) is approximately constant in the region where r(p) varies appre- 
ciably, then the near-field detection amplitude is given by 

^'„f(p„pj cxr(p,-p,). (107) 

r(p) is generally a narrow function peaked at p = 0, which guarantees that 
the positions of the down-converted photons are correlated: pj = p^. We 
note that this is a much simplified picture of near-field correlations. It has 
very recently been shown that the near-field distribution of the biphoton may 
present a complex structure due to fourth-order interference of photon pairs 
that originate from different transverse planes of the cry stal, and is strongly 



dependent upon the collinear phase mismatch [142l . |143 . 

In the far-field region, the down-converted photons exhibit ant i- correlation. 
This is understood intuitively via momentum conservation in the SPDC pro- 
cess. If the pump beam is centered around q = 0, momentum conserva- 
tion guarantees that the down-converted photons are produced such that 

^ — qfj. Since the transverse position in the far-field is associated with the 
momentum distribution at the crystal, this translates into anti-correlation in 
the detection positions in the far-field. Equivalently, if W(p) is much larger 
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than r(p), then this imphes that for the Fourier transforms of these func- 
tions the inverse is true: v{q) is much narrower than 7(g). The field operator 
corresponding to detection in the far-field as 

Eff(p) oc a f Q , (108) 



b 

so that the detection amplitude can be approximated by 

*.(P,.P.)«.(^ + ^), (109) 

where bg and 6j are scaling factors with dimension of length"^, as the ar- 
gument of the function v has units of a transverse wave vector (spatial 
frequency). Equation (11091) shows that position measurements in the far- 
field are associated with the momentum profile v{q) in the near- field. Thus, 
through position measurements in the near and far-field, one can test the 
entanglement criteria fll04p or (11051) . For example, in this case the MGVT 
criterion (11041) gives 

(A2x_)(A2p+) = A2i;A2r > 1, (110) 

where A^f and A^F are the variances of the functions v and F, respectively. 
For a Gaussian pump beam with width w at the crystal face, A^f = 



and A^F ^ ApL/(27r) 144| . where L is the crystal length and Ap is the 
wavelength of the pump beam. Since v and F are independent functions, they 
need not respect the limit on the right-hand side. In fact, for a A = 400nm 
Gaussian pump beam with width w = 1mm, and L = 1mm crystal, the 
left-hand side is on the order of 10~"^. 

The strong spatial correlation depends explicitly on the widths of the 
pump beam and phase matching function. For long crystals and strongly 
focused pump beams, this correl ation can be decreased so that the two- 



photon state is almost separable |145l . Il46| . In this regime the shape and 



form of the pump beam and phase matching function become very important. 
In the approximation where both functions are described by Gaussians, the 
two-photon state is separable when A^w = A^F. 

A condition similar to inequalities (11041) and (11051) was tested by D'Angelo 



et al. for down-converted photons using a ghost imaging setup [22|, as illus- 
trated in figure [261 Here a double slit was placed in the path of the signal 
photon, which is detected by a stationary "bucket" detector. The idler pho- 
ton is detected by scanning a point-like detector in either the near-field or 
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Figure 26: Experimental setup for test of spatial entanglement performed by D'angelo et 



al. [2^. Measurements in the near field are performed by scanning detector D2 in the 



image plane of an imaging lens system. Far-field measurements are realized by scanning 
detector D3 in the Fourier plane of the lens. Detector Di is a bucket detector, placed in 
the focal plane of a lens. 



far-field, resulting in a ghost imaging or ghost interference scenario. The 
non-classicality conditions tested were 

A^(xi — X2) < min(A^xi, A^X2); 

A^(Pi + P2) < min(AVi, AV2)- (HI) 

Violation of both inequalities by the same quantum state implies entangle- 
ment. In the experiment, A^(pi + P2) was evaluated through the visibility 
of the ghost interference setup, while A^(xi — X2) is determined through the 
ghost imaging of the double slit. By further considering the divergence of 
the down-converted beams and the localization of photons by the double slit, 
both inequalities in lllll are violated. 

7.3. Einstein-Podolsky-Rosen Non-locality 

In 1935, Einstein, Podolsky and Rosen (EPR) began the ongoing dis- 
cussion and study of quantum entanglement in their seminal paper "Can 
Quantum-Mechanical Description of Physical Reality Be Considered Com- 
plete?". EPR argued that quantum mechanics was inconsistent with seem- 
ingly reasonable notions of locality and "elements of reality" . Through analy- 
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sis of a gedanken experiment involving two particles with perfectly correlated 
position and momentum, they argued that one should be able to attribute 
well-defined values to complementary physical quantities of one of the par- 
ticles. This conflicted with the predictions of quantum mechanics, in par- 
ticular with Bohr's complementarity principle, quantified by the Heisenberg 
uncertainty relation. EPR then argued that quantum mechanics must be 
incomplete. The EPR paradox illustrates the incompatibility between what 
is now known as local realism and the standard quantum theory. 

Though EPR considered perfect position and momentum correlations, it 
is possible to apply the EPR argument to a more realistic setting. Reid and 
collaborators have shown t hat EPR-like c orrelations can be identified by the 



violation of the inequality |118l . Ill9l . Il47 



Al,^{xi\x2)AlMP2)>\. (112) 

Here A^ijj(ri|r2) is the minimum inferred variance, which represents the min- 
imum uncertainty in inferring variable ri of system 1 conditioned upon mea- 
surement of variable r2 of system 2. Explicitly, 

AL„(n|r2) = I t/6^(6)A2(ri|6), (113) 

where A^(ri|^2) is the variance of the conditional probability distribution 
V{ri\^2), which gives the probability of ri given that the measurement of 
system 2 gave result ^2, and "^(^2) is the probability that result ,^2 is obtained. 
Inequality (11121) was first violated with quadrature measurements of two 
intense beams [148]. 

Using SPDC, it is possible to produce a two-photon quantum state that 
is similar to the original EPR state. The SPDC state (see section E]) in the 
limit w — > 00 and L — > becomes 



'EPR 



j dq |qr)i l-q)^ = j \p)i \p)2 . (114) 



where the equality on the right side is obtained by Fourier transform. In- 
equality (I112p for spatial correlations was first tested experimentally by 
Howell et al. [21j, using the experimental setup illustrated in figure [271 
The near (x) and far-field (p) distributions were measured using an imag- 
ing lens system and a Fourier transform lens system, respectively. A value 
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a) 




b) 




Figure 27: Experimental setup for test of the EPR criteria performed by Howell et al. [21| . 
a) Position measurement are performed by scanning a slit aperture in the image plane of 
an imaging lens system, b) Momentum measurements are realized by scanning the slit in 
the Fourier plane of a lens. 



of A^jj^(a;i|x2)A^jj^(j9i|p2) = 0.01 << 0.25 was obtained, demonstrating the 
high correlation of the spatial degrees of freedom of the photon pairs. The 
separability condition (I104p was also tested and violated. 

A thorough discussion of the EPR paradox, incl uding the experimental re- 
sults to date can be found in a recent review article 147j . As inequality (11121) 
deals with EPR non-locality, it is generally more restrict ive than t hose which 
identify non-separability of continuous variable systems 116l IllTj . Recently, 
it has been shown that violation of the EPR inequality (I112p is sufficient to 
discard "local hidden sta te" mode l s, wh ich have been related to Schrodinger's 
"steering" phenomenon 149l . ll50L Il5l| . An EPR inequality based on an en- 
tropic uncertainty relation was introduced and tested experimentally in Ref. 
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7.^. Propagation of spatial entanglement 

The above experiments show quantum entanglement and EPR non-locality 
by identifying correlations in the near and far field of the SPDC crystal. In 
general, after free propagation, the transverse spatial correlations switch from 
a correlation in the near-field, given by equation fllOTp to an anti-correlation 
in the far-field, given by equation fllU9p . Chan et al. have shown that this 
change from correlation to anti-correlation implies that at some intermediate 
plane there rnust be very little correlation present in the spatial intensity 
distributions 144j . Following Chan et al., let us consider the case where the 
cardinal sine function in the momentum configuration is approximated with 
a Gaussian: sinc(6g^) oc exp(— afeg^), with a = 0.455(a = 0.455 optimizes 
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the approximation), and the pump is a Gaussian beam. After propagation, 
the wave function in g-space is 



4 

a„ + ib-(z) 



;il5) 



where C is a normahzation constant and a+ = w"^ / {1 + zj^^/ R"^) , a_ = aL/K, 
b+{z) = 2{z + L)/K + 2R/K{1 + zl/R^), and h^{z) = {2z + L)/K. Here 
R and z^ are the radius of curvature and Rayleigh range of the Gaussian 
pump beam in the source plane, respectively. The p-space wave-function is 
the Fourier transform of flll6p . and is given by 



— z6j 



4(a^ 



hl{z)) 
ib_{z) 



Pif 



(116) 



4(a2 +62(^)) 



One can see that both wave functions flll6p and flll7p are completely separa- 
ble (un-entangled) when a+ = a_ and hj^{z) = 



(z), and thus \E'sep( 



ipsiPsj z)ipi{pi, z). This condition depends upon the initial parameters of the 
pump beam, as well as the length of the non-linear crystal, and cannot be 
met by propagation alone. However, at a certain distance z = zq from the 



non-linear crystal, such that a+[a'^ +b'^{zo)] 
squared of the wave function is separable: 



+ b\{zo)], the modulus 



l^(Ps,P*,^o)P = \lpsips,Zo)\'^\MPr,Zo)\'^- 



;ii7) 



Thus, in the transverse plane located a distance z = zq from the crystal, there 
is no correlation present in the intensity distribution. Since at zq we have 
b+{zo){ai + b'i{zo)) ^ b^{zo){a^ + b'^{zo)), "^{p^, p^, zq) is not separable, and 
we conclude all ini tial s patial correlations have "migrated" to the phase of the 
two-photon state jl44| . Thus, to properly identify the entanglement present 
a.t z = zo, some type of interferometr ic m easurement capable of accessing 



phase information must be performed [144 



In the above arguments, the phase matching function 7(q) is approxi- 
mated by a Gaussian function, which is necessary for the separability of the 
modulus of the wave function (11171) . When considering the phase matching 
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function in the form of a sine funetion(in the momentum space), the corre- 
lations in principle never completely disappear. Nevertheless, one observes a 
strong decrease in spatial correlations at z = Zq. 

Since the correlations present in the two-photon state migrate from the 
amplitude to the phase due to local unitary evolution in the form of free 
propagation, it should be possible to retrieve the correlations by undoing the 
unitary evolution. In this vein, it is convenient to parametrize the evolution 
in a more general framework. 



7.5. Fractional Fourier Transform 

The analysis of the propagation of spatial correlations can be cast in a 
general setting using the concept of the fractional Fouri er tra nsform (FRFT) 
153l . I154I . I155I . |l56|. Since its first appearanc e in 1929 | l57l . the FRFT has 



found widespread use in qu antum mechanics |158l . Il59l . Il60l | as well as signal 
processing and optics 156|. The FRFT of order of a function i?(C) can be 
defined by the integral transform 154 



[E] {i)= jj c)E{c)dc 



;ii8) 



where the FRFT kernel is given by 



I exp [—t(f 
2tt sin d) 



exp 



I 1 h I— — - 



sm < 



(119) 



and C, and ^ are two-dimensional variables. Note that = 7r/2 corresponds 
to the usual Fourier transform. Furthermore, in the limit cf) — > 0, Fq — > 
5{C, — C)) ^-iid when cf) — )■ tt the kernel F^^ — > 5{C, + $,)■ Just as the Fourier 
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transform appears naturally in the context of Fraunhofer diffractio n, it has 



been shown that the FRFT appears in the Fresnel diffraction regime |154j . To 



identify ordinary free space propagation with the fractional Fourier transform 



it is necessary to choose properly scaled coordinates |154| . Consider the usual 
Fresnel diffraction integral corresponding to the propagation of a field from 
the plane A to plane B, as shown in Fig. [2H]a). Choosing dimensionless 
coordinates C = ^Jkiancf)/ Dp and ^ = ^Jk sin cos (j)/ Dp' , where D the 
distance between planes, leads to 

EB{i) = exp(z0) cos0exp(-«e'tan0/2)J-<^ [Ea] (^), (120) 



which expresses Eb as a FRFT of Ea multiplied by a phase term 154j |. It 



is possible to eliminate the phase term if one observes the output field on a 
spherical surface of radius = —D/ sin^ 0, as shown in Fig. [28] a). This 
spherical surface can also be mapped to a plane using a lens with focal length 
equal to — -R<^. 

It has also been shown that one can implement a FRFT using lenses 

161 . ISSj . For example, consider the symmetrical optical system shown in 



Fig. |28]b). Here / is the focal length of the lens and z is the propagation 
distance before and after the lens. Using z = 2/sin^(0/2) and / = / sine/), 
it has been shown that this optical system implements a FRFT operation 
given by 

J-^ [E] (pO = IIf^ (|p', jp) E{p)dp, (121) 

which is equivalent to the transform given in Eq. f lllSp when one chooses 
scaled variables C = ^/WfP $, = ^/kjfp'. 

A generic optical system consisting of lenses, mirrors and sections of free 
space can be described in terms of FRFTs by scaling the transverse spatial 
coordinates appropriately. The convenience of the FRFT comes from the fact 
that it can be viewed as a simple rotation in phase space. This can been seen 
most easily using geometrical optics. Consider the symmetrical lens system 
shown in Fig. [28] b). Free propagation of an optical ray r = {r,6) can be 
represented by the ABCD matrix [162] 

S^=(l 1 V (122) 



where z is the propagation distance. Passage through a thin lens is described 
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;i23) 



by the matrix 

1 

-1// 1 

where / is the focal length of the lens. Choosing z = 21 sin^(0/2) and defining 
/ = I sin as a scaled focal length, the complete optical FRFT system is given 
by the matrix 

cos (f) / sin I 
— J sin cos 4 



SzLiSz 



(124) 



Matrix fll24p represents a 0-order FRFT, and is recognized as a rotation 
matrix scaled by /. This scaling is necessary since r has dimension of length 
and 6 is adimensional, and thus the scaled focal length / acts on 6 so that 
f6 has dimension of length. If the scaled focal lengths / of two FRFTs are 
equal, it is easy to see from the ray matrix in Eq. f ll24p that the FRFT 
is additive, that is F^^F^^ = F^^+^j, such that the order of the combined 
optical FRFT is 0i + 02- By introducing the FRFT with properly scaled 
coordinates, and using the additivity property, one can describe propagation 
through any first-order optical system consisting of lenses, mirrors and free 
space. Propagation is now parametrized solely by the order of the overall 
FRFT, given by the angle 0. 

In q uantum-mechanical operator formalism, the FRFT operator is defined 



as 



156 



e^*/2 exp 



(125) 



where x and p are the dimensionless position and momentum operators sat- 
isfying [x, p] = i. This operator is equivalent to the evolution operator of the 
quantum harmonic oscillator, which has the hamiltonian H = (x^ + p^)/2. 
Application of the FRFT operator to the x and p corresponds to rotation 



of angle in phase space |16ll . |156 



cos (p sm < 
— sin 6 cos I 



(126) 



where x^ and p^ are the rotated operators. 

The FRFT is a rotation in the x,p phase space. In this respect, it can be 
used to measure the marginal probability distribution along any ax is in phas e 
space, and to perform tomography of the spatial Wigner function 1631 . Il64 . 
We note that spatial tomography of a field can also be perfor med using 
interference of rotated and displaced halves of the field 



be p 

m 



16a 166 . 
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1.6. Correlations at Intermediate Planes 

To describe the spatial correlations of photon pairs at some propagation 
distances Zg and Zi, through arbitrary first-order optical systems, it suffices 
to consider FRFT's of arbitrary order angle a and /3. A two-photon state 

The 



:(2) 

/3 



after arbitrary propagation is given by \'^a,i3) = 
two-photon wave function then becomes '^a,i3{xi,X2) = {xi,X2\'^a,i3), where 

^Q,/3(a;i,X2) = jj dx[dx'2 {xi \ Fa \x'i) {x2\ Ixg) ^(x^jXa) (127) 

= jj dx[dX2Fa{Xi,x[)Fp{x2,X2)'^{x[,X2), 

and the kernels are defined in Eq. flllQp . For simplicity, we consider only 
one spatial dimension. As an example, let us consider the propagation of the 
EPR state in one dimension 



EPR 



> 



dxidx2S{xi — X2) \xi)^ \X2) 



2 ' 



(128) 



which presents a perfect correlation, since detection of photon 2 at position 
X projects photon 1 onto a position eigenstate \x). This situation is approxi- 
mated by the state produced by SPDC whe n the pump beam can be treated 
plane wave. This state evolves to |167l | 



\^^^) = Aa+p jj dxi(ix2exp 



.cot (a -|- (3) 



{pI + pD 



exp 



—I 



Xi ■ X2 

sin(a -|- (3) 



(129) 



We can now analyze the type of correlations present. Whenever a + 
13 = (mod27r), the original state fll28p is recovered. That is, the EPR 
state (11281) is an eigenstate of operators of the type fa^2n-a, ^a^An-a, etc. 
When a + (3 = 71 (mod27r), the correlated EPR state (11281) evolves to an 
anticorrelated EPR state 



dxidx2S{xi + X2) \xi)^ 1x2)5 



(130) 



In this case the detection of photon 2 at x projects photon 1 onto the state 
|— x). When a + (3 = n/2 (mod27r), this state becomes 



\Q) 



dx |x)^ \p{x)) 



2 ' 



(131) 
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which presents no intensity correlation. Here \p{x)) oc J exp{ip{x) ■ x) \x) is 
the momentum eigenstate conjugate to \x). An equivalent result is found 
for a + P = Sir/ 2 (mod27r). We note that the conditions for correlation, 
anti-correlation, and no-correlation depend on the sum of the FRFT angles 
of the down-converted fields, and not the individual angles a and (3. Thus, 
for any propagation characterized by an FRFT J-'q on photon 1, one can 
find a suitable transformation J-"^ on photon 2 such that a correlation, anti- 
correlation or no intensity correlation is recovered. These conditions were 



tested experimentally in Ref. [167 



This simple picture drawn for the ideal EPR-state is followed approxi- 
mately by the two-photon state in Eq. (11171) . Let us consider two simple 
cases. First, let us assume that the FRFT angles are the same for both 
down-converted photons: a = p. Then, no intensity correlation is observed 
whenever 

2 a^{al + bl) -a+{al + bl) 

cot a = , (132) 

a+ — a_ 

where parameters a± and b± are defined in flll6p and b± refers to b±{z) at 
the crystal plane. Consider another simple case in which 6+ = 6_ = at the 
crystal face, with no restriction on a and f3. Analytical calculation shows 
that, in order to have no intensity correlation, 

cot a cot /3 = a_a+. (133) 

Eq. (11331) is satisfied by FRFT orders such that a + f3 = it/ 2 (mod27r) or 
a + f3 = 37r/2 (mod27r) only when a_ = l/a+. Nevertheless, the intensity 
correlations present in the state Eq. (I117P propagate in a fashion similar to 
idealized case of the EPR state. 

The above arguments apply to intensity correlations, the absence of which 
does not necessarily imply that a quantum state is separable. Let us now 
discuss the identification of entanglement via intensity correlations of trans- 
verse spatial variables. The DGCZ inequa lity (11051) for transverse variables 



rotated using the FRFT operator (11261) is 168 



((AX'_)2) + {{AP'f) = i±££^f^±« [((AX_)2) + ((AP+)2)] + 



2 

lz£2E^ K(AX+)2) + ((AP_)2)] - 
5ii^[({X+,P+})-2(X+)(P+)] + 
2H^[({X_,P_})-2(X_)(P_)], (134) 
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Figure 29: Laguerre-Gaussian and Hermite-Gaussian modes. 



where X'_ = Xq, — and P'_,_ = p^, + p^, are defined in terms of tlie rotated 
variables, and X_ = Xi — X2 and P+ = Pi + P2 in terms of tlie variables at 
the source. The sum of variances for the rotated variables coincides with the 
sum of variances for the variables at the source when a + /3(mod27r) = 0. 
If the two-photon state is entangled, with intensity correlation in the x and 
p variables at some initial plane, then, for any propagation of the signal 
photon, characterized by a, it is possible to find a propagation of the idler 
field, characterized by /3, so that an intensity correlation is recovered and 
entanglement can be identified. It is important to note that Eq. (11341) does 
not depend on the state and is applicable to any bipartite continuous variable 
systems. Eq. fll34p was experim ental ly tested for intermediate propagation 



planes of SPDC photons in Ref. [168 



8. Transverse Modes 

A simple way to describe spatial correlations in the paraxial regime is 
through the use of transverse spatial modes described by discrete indices, such 
as the Laguerre Gaussian (LG) modes or the Hermite Gaussian (HG) modes. 
Like the well known Gaussian beam, the Hermite-Gaussian and Laguerre- 
Gaussian beams are also solutions to the paraxial Helmholtz equation (HTi) 



35l |. In section IHTTI we will introduce the Laguerre Gaussian modes, which 
have well defined orbital angular momentum (0AM). In section [8^ we review 
the conservation of 0AM in SPDC and the generation of entangled 0AM 
states. In sections 18.51 and 18.61 we present the Hermite-Gaussian modes and 
their application to SPDC. 
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8.1. Laguerre Gaussian modes: orbital angular momentum of light 

Since the 1990's it has been known that any electromagnetic paraxial 
beam with an azimuthal phase depen dence of the form e^^^ carries an or- 
bital angular momentum Ih per photon jl69|. The Laguerre- Gaussian beams 
are probably the most well-known and studied examples of beams carrying 
orbital angular momentum. They are given by 170 



'w{z) \ w{z) j ^ \w{zy 



p 



exp < —t 



2R 



{n + m + l)'y{z) 



w{zy 



(135) 



where (p, 0, z) are the usual cylindrical coordinates, is a constant which 
depends on the azimuthal index i and radial index p. are the Laguerre 
polynomials. Here z is the longitudinal propagation direction, R{z) is the 
radius of curvature, w{z) is the beam waist, and 7(z) is the phase retardation 
or Gouy phase. The parameter z^ is the Rayleigh range. The order of the LG 
beam is defined as A/" = |£| +2p. The usual Gaussian beam is the zeroth-order 
Uq beam. See plots of LG modes in Fig. I29r ieft). 

In addition to interesting implications in classical and quantum optics, 
the orbital angular momentum of a light field raises possibilities for techni- 
cal applications. It has been shown that the orbital angular m ome ntum 
of light can be used to rotate micro-particles in optical traps |17l| . In 
terms of quantum optical applications, the orbital angular momentum of 
single photo ns in LG mod e s provide s a p ossible d- dimensional qudit encod- 



ing scheme [172, Il73l . 1 174 . Il75l . Il76l . Il77l |. which allows for the creation of 



multi-dimensional entanglement in discrete bases. Devices that discriminate 
the orbital angular momentum of Lag uerre-Gau ssian beams have been pro- 
posed |l78| and experimentally tested 



179, 180 



.2. Entanglement and Gonservation of Orbital Angular Momentum in SPDG 

The first experimental investigation of 0AM conservation in SPDG sug- 

in contrast to previous experimental 



gested that QAM is not conserved [181 



investigations of other non-linear optical processes such as second-harmonic 
generation (SHG). Using an LG mode as the fundamental beam in SHG 
showed a conservation relation, such that the QAM of the second-harmonic 
beam is double that of the fundamental beam: 



-SH 



21 



fund 



182, 183 . In 
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these experiments, the intensity and phase structure of the output beam was 
observed using a CCD camera. The same measurement method was used 
in the inve rse e xperiment, investigation of 0AM conservation in SPDC by 
Arlt et al. 
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However, by imaging the near and far field of the coUinear 
down-converted fields using a CCD camera, correlations in second-order are 
observed. Theoretical investigations by Arnaut and Barbosa showed that, 
under suitab le ex perimental conditions, 0AM conservation appears only in 
fourth-order 172| . The two-photon state has well-defined 0AM equal to that 
of the pump beam, and 4ignai+^idier = -^pump- Moreover, they showed that the 
two-photon state is entangled in 0AM, and thus the individual signal and 
idler beams do not have well defined 0AM. For this reason, the experiment 
by Arlt et al. did not observe 0AM conservation. Shortly thereafter, a clever 
experiment by Mair et al. confirmed the conservation and entanglement of 
0AM in SPDC [l8^ . 

The experiment by Mair et al. 



184| is illustrated in figure 



The 

0AM correlation was measured using forked holographic masks and single- 
mode optical fibers to project onto individual 0AM values of m. The forked 
hologram mask functions in much the same way as a diffraction grating. A 
Gaussian mode passing through a mask with m dislocations and diffracting 
into the n^'^-order becomes an LG mode with / = mn. The detection scheme 
of Mair et al. exploits this process in the reverse. A mode with 0AM m 
that diffracts into the first-order of a forked mask with m dislocations will 
now be in a Gaussian mode. Single mode fibers are used to select only those 
photons in the Gaussian mode, which are then registered by single photon 
detectors. The pump laser beam with / = —1,0,1 was shown to produce 
entangled photons with rrisignai = l — ^idier- However, this confirms only a 
classical correlation between the 0AM of the down-converted photons. To 
prove that the photons are indeed entangled in 0AM, Mair et al. projected 
the photons onto superpositions of 0AM states. This was done by shifting 
the position of the holographic mask. A shifted mask with m dislocations and 
the single-mode fiber detection scheme described above projects the single- 
photons onto superposition states of the form a |0) + /3 |m), where a and (3 
depend on the position of the mask. 

The two-photon state for frequency degenerate down-converted fields 
(As = Aj = 2Ap) in a quasi-coUinear ge ometry can be dec omposed in terms of 



LG modes in a relatively simple form jl72l Il85l . Il86l . Il87j . The wave function 
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Figure 30: First measurement of 0AM entanglement in SPDC by Mair et al. ^184,]. Forked 
holographic masks of 0AM Ig and li are used in conjunction with lenses and single mode- 
optical fibers to project onto modes with well-defined 0AM. To project onto superpositions 
of 0AM, the masks are shifted in the transverse direction. 



in this case can be written as 



^{Ps, Pi) = U„ I — I -A 



V2 



V2 



(136) 



where Up describes the pump beam profile and J-" is the phase- matching 
function. 

Decomposing this wave function into signal and idler LG modes results 
in an alternative form: 



^iPs,p.) = Y.Y.K:^K:(ps)K\(p^ 



(137) 



where the coefficient A^i.^'l! is given by |186 



PsPi 



rl I Ps + Pt 
's^'t^i ^p 



V2 



7 



(138) 
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Assuming that the crystal length L is much smaller than the Rayleigh range 
zr of the pump beam, and the function ~ 1, the coefficient N^^j^, becomes 

Ki. « ^^s+ia [qdq <(v^g)<=(g)<'(g), (139) 



where Vp{q) is the Fourier transform of the radial component of U!p{p). The 
presence of the Kronecker delta function guarantees that 0AM is conserved: 
k + ls = l-i and the down-converted signal and idler photons are emitted with 
correlated 0AM. For simplicity, let \m) denote a single photon state carrying 
an orbital angular momentum mh and I be the azimuthal index of the LG 
pump beam. Coherence of the SPDC process guarantees that t he tw o-photon 



state can be written as a Schmidt decomposition of the form |145 



|^)= Pn.\l-m)\m). (140) 

m=— oo 

Here the states |m) represent spatial modes with well-defined 0AM m. 

The derivation of the two-photon state (I140p is valid for quasi-collinear 
paraxial fields in the thin-crystal approximation. Other effects, such as the 
birefringence and anisotropy of the non-linear medium, as well as a non- 
collinear geometry, were ignored. Some authors have studied SPDC and 



0AM conservation in more general contexts [45|, |43|, |188| . A c omprehensive 
discussion has been provided recently by Osorio et al. [l88j. They have 
shown that when one considers all the photon pairs produced by the crystal, 
distributed along the entire down-conversion cone, 0AM is conserved. In 
a more realistic setting, in which signal and idler photons are detected in 
small regions of the down-conversion cone, the two-photon wave function 
may be spatially asymmetric and not reproduce the angular spectrum of 
the pump beam. The deformation generally corresponds to an ellipticity 
that depends upon the non-collinear angle y9, the crystal length, and the 
width of the pump beam w , and can be quantified by the quantity = 
w/siinp. When the crystal length L <^ L^, the two-photon wave function 
reproduces the transverse field of the pump beam, and 0AM is conserved. 
If this condition is not fulfilled, the two-photon wave function is not given 

by ^(p^,Pj) = Up and 0AM will not be conserved. Ellipticity of 

the two-photon spatial wave function caused by non-collinear geometry and 



focusing of the pump beam has been observed experimentally |47|]. Torres, 
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Osorio and Torner 189| have considered the extreme situation in which the 
down-converted fields are either counter-propagating and perpendicular to 
the pump field, or co-propagating with respect to each other but counter- 
propagating to the pump field. In these cases the usual 0AM selection rule 
(as in the collinear geometry above) does not apply. 

Another consideration is the spatial walk-off caused by the birefringence 
of the non-linear SPDC crystal, which occurs for either the pump and/or one 
of the down-converted fields. In the case of type-I down-conversion, in which 
only the pump field is polarized in the extraordinary direct ion, the critical 



parameter is the "walk-off length" L^, = w / tanpo |43|, Il88l |. Here po is the 
spatial walk-off angle, which depends upon the extraordinary refractive index 
and the propagation direction in the crystal. When the length of the non- 
linear crystal L <C L^, the effect of walk-off on the shape of the two-photon 
wave function is negligible, and 0AM is conserved. Fedorov 4^ . 



8.3. Probing the phase structure of the two photon state: two-photon inter- 
ference 

Two-photon interference at a be am splitter was first demonstrated by 
Hong, Ou and Mandel (HOM) jl90 |. It has subsequently become the key 
process for two-photon quantum logic, and is u s ed in quantum information 
protocols such as B ell-state measurements 19ll . Il92[ | and to construct two- 



photon logic gates |193l . Il94j . Il95l . Il96l . Il97l |. On the one hand, these ap 



plications are designed for perfect interference with a single spatial mode. 
On the other hand, multimode interference is a use ful metho d to probe the 
transverse phase structure of the two-photon state [l98|, [l99|, and has been 



used to show the conservation of orbital an gular momentum in SPDC [186 



as well as to quantify spatial entanglement 20d |. It has also been shown to 
be useful in Bell-state analysis of polarization entangled photons 20 1| . and 
to produce two-photon fields exhibiting spatial antibunching 113 . 



Let us illustrate the probing of the phase structure. Consider the HOM 
interferometer shown in Fig. [311 in which two photons are generated by 
SPDC and then refiected onto opposite sides of a beam splitter (BS). If 
the path length difference is greater than the coherence length of the down- 
converted photons, then there is no interference and the photons leave either 
side of the beam splitter randomly. Here we will assume that lengths of 
paths s and i are equal. Using the reference frames illustrated in Fig. [311 
the annihilation operators in exit modes 1 and 2 after the beam splitter can 
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Figure 31: Illustration of Hong-Ou-Mandel interference with multiple spatial modes. 



be expressed in terms of the operators in input modes s and i: 

ai(q, a) = tasiqx, Qy, o") + irdiiq^, ^r) (141a) 
a2(q, cr) = tai{qx, qy, a) + iras{qx, -qy, o"), (141b) 

where a is the polarization, qx and qy are the transverse components of 
the momentum and t and r are respectively the transmission and reflection 
coefficients of the beam splitter (|tp + |rp = 1). A field reflected from the 
beam splitter undergoes a spatial inversion in the horizontal {y) direction, 
while a transmitted field does not suffer any inversion, as illustrated in Fig. 
[3n The negative sign that appears in the qy components of the reflected 
field is due to this spatial inversion. The two-photon wave function is split 
into four components, accordi ng to the four possibilities of transmission and 
reflection of the two photons |198l . [l99|: 

(r+ , r _ ) = (r+ , r_ ) =itr exp j {x'i + yl 



Z 

[W {x+, Z) + W {x+, Z)] , (142) 
*i,(r+,r_) = t^exp |^(x2 + y'_)]w {x^,y+, Z) , (143) 
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*„(r+, r_) = -r^ exp { '^{xl_ + yl) \W (a;+, -y+, Z) , (144) 



where r± = {x±,y±) = {xs/2 ± Xi/2,ys/2 ± yi/2). Let us assume that 
t = r = I/V2. If the pump beam profile W is symmetric with respect to 
the y coordinate (W[|/] = W[— y]), then the two photons always leave in 
the same output port of the BS, since the amplitude for coincidence counts 
in different output ports is ^tt + ^rr = 0. This results in the usual HOM 
interference dip. If W is antisymmetric with respect to y {yV[y] = — VV[— y]), 
then "^tr = ^rt = 0, and the photons always leave the BS in different output 
ports, giving an interference peak. 

The analysis above considered that all other degrees of freedom are in 
symmetric states. Combining the spatial degrees of freedom with the polar- 
ization degree of freedom, for instance, provides a way to control two-photon 
interference of different polarization states. This allows for partial Bell-state 
analysis in the coincidence basis \20l\ . as well as for the creation of a two- 
photon "singlet" beam, in which the photon pairs in the same beam are in 
an anti-symmetric polarization state, and consequently are spatially anti- 



bunched 113 



Other experiments have used two-photon interference to probe the phase 
structure of the two-photon state, an d have sh own that the down-converted 



photon pairs are entangled in 0AM |186l Il99| . For a pump beam in a LG 



mode, the two-photon coincidence detection probability at the output of the 
HOM interferometer is 

P{r,,r2)<x\ul{R)\^ single, (145) 

where Up{R) is the radial component of the LG mode, i is the 0AM of the 
pump beam which is transferred to the two-photon state, R = |Pi + P2I) 
and sin6' = (pi sin^i + p2 sin 02) /-R- When r2 = 0, the coincidence profile is 
P(ri,0) oc |Up(pi/v^)psin^£0i. Thus, fixing one detector and scanning the 
other, os cillat ions with period i will appear in the azimuthal coordinate 0i. 
In Ref. fl86| . the phase structure was observed for LG pump beams with 



£=1,2. 

8.4- Quantifying Spatial Entanglement 

In section [71 the detection of spatial entanglement using entanglement 
witnesses was discussed. In many cases, these witnesses involve only a few 
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measurements of the transverse spatial variables. Quantifying spatial entan- 
glement, on the other hand, presents a considerable challenge, due to the 
large dimensionality of the involved Hilbert spaces. Nonetheless, theoret- 
ical and experimental methods exist which allow for the quantification of 
bipartite spatial entanglement of photon pairs. 

8.4-1 ■ Schmidt decomposition of SPD C 

The Schmidt decomposition 202| of the t wo-ph oton state produced by 



SPDC has been derived by Law and Eberly [145|. Absorbing any phase 
factors into the s tates \m) allows one to write the two-photon amplitude 
(11401) in the form [145 



(146) 



-oo n=0 



where the real numbers Ar 
tors. The functions 
defined 0AM 



145 



v„ m are the eigenvalues of the reduced density opera- 
are normalized transverse mode functions with well- 
The entanglement depends upon the number of non-zero 
Schmidt coefficients \n,m,- Each \n,m depends on the 0AM wavenumber 
and width of the pump laser beam, and the length of the non-linear crystal, 
and can be calculated by determining the overlap integral of the product of 
mode functions with the two-photon amplitude given in equation ( 1136^ : 



PsPiU 



nl- 



[PsK. 



V2 



7 



(147) 



Again, effects due to crystal anisotropy and non-collinear geometry have 
been ignored. Law and Eberly 145| have used the decomposition (I146P to 
calculate the entanglement of the two-photon state generated by a Ga ussian 
pump beam using the so-called Schmidt number, 5* = |l45 |. 
They have shown that the entanglement depends upon the parameter L/2zr, 
where zr is the Rayleigh range of the pump beam and L is the crystal length. 
If one approximates the phase matching fun ction by a Gaussian function with 
the same variance, 7(g) ^ exp(-LgV8ir) [lii], the Schmidt number in this 
case can be put in the analytical form S = {^J L/Azr + a/4z/j/L)^/4. The 
entanglement in this Gaussian approximation serves as a lower bound, since 
Gaussian states are alwa ys le ss entangled than non-Gaussian states with the 
same covariance matrix 203 . 
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For typical values of the experimental parameters, the predicted Schmidt 
number can be very large. For example, for L = 1mm and Zji = Im, the 
predicted Schmidt number is on t he o rder of 10^. However , as has been 
pointed out by van Exter et al. 204j . and Osorio et al. 188l |. even in 



the case of coUinear propagation and negligible crystal anisotropy, this does 
not necessarily correspond to the usable entanglement present in a realistic 
experimental scenario. The effect of spatial filtering due to the finite spatial 
bandwidth of an optical setup is to decrease the Schmidt number, and hence 
the entanglement. This can be a rather large effect, of several orders of 
magnitude. 

In Ref. [205^, Di Lorenzo Pires et al. demostrated that for the case of a 
two-photon pure state such as (jSl]), the overall Schmidt number (including 
contributions from radial and azimuthal parameters) is given by the inverse 
of the overall degree of coherence of either down-converted field. Assuming 
that-say-the signal field is quasihomogeneous, the Schmidt number takes the 
form 

where and Iff are the near- field and far-field intensity distributions of 
the signal photon field, respectively. Thus, through second-order intensity 
measurements, the spatial entanglement of the two-photon state can be de- 
termined. It was shown that the entanglement depends critically on the 



collinear phase mismatch |205| . Schmidt numbers above 400 were obtained 
experimentally. 

An experimental technique to measure the Schmidt number associated 



with the 0AM decomposition (I146p was developed by Peeters et al. [200 
who investigated multimode HOM interference in the presence of image ro- 
tation of the signal beam before the beam splitter. They showed that the 
amount of spatial entanglement can be inferred from the visibility of the 
HOM dip: 

oo 

¥{;&)= P^cos(2mt9), (149) 

m=— oo 

where d is the angle of image rotation of the signal field. Here Pm are the 
coefficients of the Schmidt decomposition f ll40p . To quantify spatial entangle- 

oo 

ment, an azimuthal Schmidt number can be defined as Saz = ( Yl ^m)^^- 

•m=—oo 

The Schmidt coefficients can be determined by measuring V{d) for a number 
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of values of i!} and inverting the Fourier series fll49p . In Ref. [200| . this was 
done for a number of different types of spatial filtering systems. Schmidt 
numbers as high as 7.3 ±0.3 were determined experimentally. More recently, 
Di Lorenzo Pires et al. performed a sirnilar experiment in a coUinear ge- 
ometry using a bucket detector system 206| . This impressive experiment 



allowed for the measurement of the complete 0AM distribution of the cor- 
related down-converted photons, and Schmidt numbers above 35 were mea- 
sured. This was in excellent agreement with theoretical predictions based on 
the experimental parameters. It was also shown that the phase mismatch 
could be used to increase the entanglement. In this case, the mismatch could 
be controlled through the temperature of the periodically poled crystal. 

Another method to quanti fy en tanglement using multimode HOM inter- 
ference has been proposed in [l46|. In this scheme, it is possible to meas ure 



the concurrence directly using two copies of the quantum state |207l |208 



An interesting series of experiments have been performed employing phase 
late s arid sin gle-mode fibers to detect and analyze spatial entanglement 
209l . I210I . l211j . Two main types of phase plates have been investigated: spi- 



ral phase plates and Heaviside step plates. When a Gaussian mode passes 
through a half-integer spiral pha se plate, it is transformed into a superposi- 
tion of LG modes, given by |212| 



(150) 



where fp/{a) is a function which depends upon the orientation angle a of the 
phase plate. By considering the inverse scenario, it is apparent that the spiral 
phase plate transforms the superposition (11501) into a Gaussian mode, which 
can be selected with a single mode optical fiber. By changing the orientation 
angle a of the plate, one can project onto a number of superposition states of 
the form fll50p . In Ref. 212| it was shown that the number of non-zero terms 
in the superposition fll50p can be quite large (~ 10^ — 10^), and depends upon 
the half-integer value of the plate. 

Half-integer spiral phase plates were used to observe fractional 0AM en- 



tanglement in down-converted photons |209l |. The experimental setups are 
similar to that of Figure |30l with the forked holographic masks r eplac ed by 



the phase plates. The coincidence detection probability, given by |209 



a,- oc vr 



Oii\y 



(151) 
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depends only upon the relative angle \as — between the analyzers in the 
signal and idler fields, which is typical of a quantum correlation. 

Another type of phase plate which has been studied is a Heaviside phase 
plate, in which there is relative pha se di f feren ce of vr between the phase im- 
parted by each section of the plate |212 . 210 1. The transmission function is 



given by T{(f), a) = 1 — 2[0(0 — a) — 0(0 — a — /?)], where 6 is the Heaviside 
step function, a is the orientation angle of the plate and /3 is the angular 
width of the vr-step section. The simplest case corresponds to two half-circle 
sections (/3 = tt). These step phase plates can also be used to probe multidi- 
mensional entanglement. Pors et al. have shown that the amount of usable 
entanglement that i s availabl e is given by the dimensionality of the mea- 
surement apparatus 210l . |211| . The authors have dubbed this the "Shannon 



dimensionality" , i n an alogy with the Shannon number of a classical commu- 



nication channel |213l |. The Shannon dimensionality is an effective Schmidt 
number, corresponding to the amount of entanglement which can be accessed 
by the measurement system. In the case of the phase plate analyzers, the 
dimensionality has a simple experimental interpretation. It is given by 27r 
divided by the area under the normalized coincidence curve, plotted as a 
function of the relative angle of the phase plates. Experimental and theoret- 
ical results show that the half-circle plate corresponds to dimensionality of 
three, while a quarter section plate gives dimensionality of six. The authors 
also considered step plates with multiple sections, and have shown that it 
should be possible to probe spaces with dimension on the order of about 50 



211 



There is thus an inconsistency between the amount of entanglement which 
is theoretically available and the amount of entanglement actually accessible 
in a given experimental scenario, due in particular to the geometry of the 
down-conversion setup and the spatial bandwidth of the optical systems used 
in manipulation and detection of the fields involved. 

Another experiment investigating this kind of entanglement quantifica- 
tion, used entangled states prepared placing apertures in the path of the 
beams generated by SPDC. In t his case the aperture photon path defines 
the qubit state. In reference |214| . two double-slits placed at the path of the 
idler and signal photons were used to generate two-qubit entangled states. 
States with different degrees of entanglement were prepared by modifying the 
pump beam profile that generates the twin photons. Measurements of either 
two-photon conditional interference or marginal probabilities were used for 
characterizing entanglement in two-qubit spatial pure states. 
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The concurrence C is an entanglement quantifier for bipartite systems 



215| . For a general pure entangled state of two qubits, 

= cnlll) + ciollO) + coilOl) + CoolOO) , (152) 

the concurrence is given by 

C(?/') = 2|ciiCoo - ciocoil . (153) 

State fll52p may be written as a Schmidt decomposition |216 

|^)=c;jll)' + c'oo|00)', (154) 

where {i = 0, 1) is the Schmidt basis, and the coefficients c'^^ are real 
and positive. In terms of these coefficients, the concurrence is given by the 
simpler expression 

Ci^P) = 2c[,c',, . (155) 



In Ref. |214| . states |0), |1) represent the single photons which pass through 
either slit of the double slit aperture. The strategy used for characteriz- 
ing entanglement consists in performing measurements in the Schmidt basis, 
either by directly measuring the Schmidt coefficients or through a direct 
connection between concurrence, given by Eq. ( 11551) . and the visibility of 
interference patterns obtained when one of the detectors is used as trigger. 
The Schmidt basis is selected by letting the photons propagate through an 
appropriate lens system and keeping one of the detectors fixed at some spe- 
cific positions, while the other is scanned at the Fourier plane of the lens. 
The Schmidt coefficients are then obtained from these conditional interfer- 
ence patterns. In the second method the fixed detector used as trigger is 
completely opened detecting all photons that crossed one of the double-slits 
while the second point detector is scanned at the lens Fourier plane. This 
corresponds to a measurement of the marginal probability which is defined 
as the probability, Pi{xi), of observing a photon idler at Xj and the signal at 
any location (— oo < < +oo) js^] and is related to the concurrence as 



Pi(x,) oc 1 + VI - [C(*)]2cos(aXi) . (156) 

In reference [55] different two-qubit states are generated by using collinear 
type II SPDC and a double-slit. Different states are produced by tailoring 
the imaging system in between the two-photon source and the double slit. 
More than 30 different two-photon states are characterized by measuring 
their complete two-photon interference patterns in the far field of the double 
slit. 
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8.5. Hermite- Gaussian modes 

The mode decomposition of the two-photon state produced in SPDC can 
also be performed using the Hermite-Gaussian modes, which are given by 
the complex field amplitude |170 



Unmix, y, Z) —C„ 



exp 



w{z) \w{z) 
2R{z) 



V2y_ 

w{z) 



exp 



— (n + m + 1)7(2:) 



(157) 



Here Cnm are the coefficients, Hn{x) is the n*'^-order Hermite polynomial, 
which is an even or odd function of x if n is even or odd, respectively. The 
order M of the beam is the sum of the indices: M = m + n. The usual 
Gaussian beam is the zeroth-order beam. See plots of HG modes in 
Fig. E^right). 

8.6. Down- conversion with Hermite-Gaussian modes 

Consider SPDC in which the pump beam is des c ribe d by a HG mode 
U^^. The two-photon state can be expanded as |l46l . 1217 



E 

j,k,u,t=0 



/~i nm 

^jkut 



Vjk) \Vut) 



(158) 



where \vab) are Hilbert space vectors such that Vab{q) = {Q\vab) corresponds 
to the HG mode written in fc-vector space, given by the Fourier transform of 
Eq. f ll57p . The coefficients Cp,I^f are related to the pump beam and phase 
matching function through 

where ^nm{<is^ Qi) is given by the product of the angular spectrum for a HG 
pump beam and the phase matching function (see section [3]). An analytical 



expression for the coefficients CJI^f^ was calculated in Ref. [217J, showing 
that spatial parity is conserved in the x and y directions. That is, the indices 
j + u {k + t) must have the same parity as the pump beam index n (m). Since 
the parity of each down-converted photon is undetermined, yet the sum of 
parities is well-defined, the two-photon state is entangled in parity. Thus, 
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assuming that the pump beam has well-defined parity in the y direction, and 
ignoring the x direction, the two-photon state can be written as 



' even 



= ^(|V'ee) + |V'oo)) 



or 



/odd 



V2 



(I^EO) + I^Oe)) 



(160) 



(161) 



where E and O stand for the y indices of photons 1 or 2 being even (odd) 
and odd (even) numbers, respectively. If the pump beam is a HG mode, the 
states IV'xy) (-^j ^ = E,0) are given by linear combinations of modes with 
the same y parity. For example. 



(162) 



k even, t odd j,u 



and similarly for the other states in fll60p and f ll6ip . Yarnall et al. 2181 . 1219 
have used this type of parity entanglement to violate Bell's inequahties with 
parity measurements, which we describe in detail below in section 18.71 

Using the expansion f ll58p . the concurrence of the two-photon state can 
be calculated as a function of the pump beam profile and the length of the 
non-linear crystal. Numerical results reported in |l46 | show that in general 
the entanglement increases with the order n -|- m of the pump beam. 

8. 7. Bell Non-locality 

The Einstein-Podolsky-Rosen paradox showed an inconsistency between 
local-realism and the completeness of quantum mechanics. In 1965, John 
S. Bell con sider ed the predictions of local realism against those of quantum 
mechanics 220|. His argument, made in the context of David Bohm's version 
of the EPR paradox involving two entangled spin-1/2 systems 22 1| . showed 



a conflict between the predictions of quantum mechanics and local realistic 
theories. This provided an experimentally accessible testing ground for these 



theories. The Clause r-Ho rne-Shimony-Holt inequality |222j . also derived in 



dependently by Bell 223| . shows that any local realistic theory should obey 



S = \E{a, h) + E{a', h) + E{a, h') - E{a', h') \ < 2, 



(163) 



where —1 < E{a,b) < 1 is the correlation function of dichotomic measure- 
ments a and b on systems A and B. For a pair of maximally-entangled spin-^ 
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particles, quantum mechanics predicts S'max = 2-\/2. This set the stage for 
tests of quantum mech anics versus local realism beginning with the work of 
Freedman and Clauser 224| and Aspect and collaborators 2251 . |226| . With 
the advent of down-conversion sources, recent experiments have reached a 
considerable degree of sophisti cation, wi th the objective of performing a con- 
clusive loop-hole free test [48|, I227L l228| , and also to discard more complex 



alternatives to standard quantum theory |229|, l230 



An early adaptation of the Bell inequality to continuous variables was 
provided by John Bell, who considered linear combinations of position and 
momentum observables of the form x-|-ap, where a is a continuous parameter 



223[. By binning the results, for instance into non- negative and negative 
values, one can translate these continuous measurements into two outcomes 
and apply the Bell-CHSH inequality fll63p . Bell reached a striking conclusion: 
quantum states with positive Wigner functions should allow for local realistic 
descriptions, since in this case the Wigner functi on it self serves as a classical 



probability distribution of the x and p variables [223[. A sur prising example 



is the EPR state (11281) . with Wigner function given by 223 



Wepk{xi,Pi,X2,P2) = 27l6{Xi - X2)S{pi +P2) 



(164) 



which is non-negative everywhere. This state displays an infinite amount of 
entanglement, and as such it is remarkable that it does not violate the Bell- 
CHSH inequality (I163p . The EPR state is not physical, however it represents 

the limit a_/aH > (and 6+ = 6_ = 0) of the Gaussian state (11171) . 

Despite these arguments against the use of continuous v ariables in this 



context, it was later shown by Banaszek and Wodkiewicz |23lL |232| that 
the EPR state does violate the Bell-CHSH inequality when one considers 
measurements of the displaced parity operator. The parity operator does 
not have a well-behaved Wigner- Weyl representation, and thus expectation 
values of this operator do not correspond to averages over classical probability 
distr ibutions. This point has been explored in more detail by Revzen et al. 
233| . Thus, the state of a given system alone is not enough to determine 
whether it is suitable or not for a Bell-type experiment, but rather one must 
also take into account the observables measured. 

Extending the idea of the parity operator, it possible to define a set of 
pseudo-spin operators, forming an analogy with the Pauli operators for a 
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Figure 32: Asymmetric interferometer used as a parity projector. When 61 = 62 and the 
beam splitters are 50:50, transverse modes with even parity leave the interferometer in 
port 1, and odd transverse modes in port 2. 



spin-1/2 system. Chen et al. have shown that the operators 234 



X = ^ |2r2 + 1) {2n\ + \2n) {2n + l\ 

n=0 

00 

Y = 2 ^ |2n + 1) {2n\ - \2n) {2n + 1\ 

n=0 

00 

Z = Y^ \2n) {2n\ - \2n + 1) (2n + 1| , 



(165a) 
(165b) 
(165c) 



n=0 



defined in terms of photon number Fock states \n), obey a SU(2) algebra. 
One can recognize Z as the parity operator. These operators were originally 
formulated for single-modes of the electromagnetic field, and their realization 



presents considerable experimental challenges |234 



An insightful contribution was made by Abouraddy et al., who recog- 
nized that the opera tors (11651) could be constructed for the spatial degrees of 
freedom of photons 218| . This is possible due to the isomorphism between 
the photon number Fock states |n), and the spatial degrees of freedom of a 
field described by a one-dimensional Hermite Gaussian mode with index n 
218| . To employ the full power of the analogy between the spatial version 



of the pseudo-operators (11651) and the Pauli operators, one must be able to 
implement arbitrary rotations, composed of linear combinations of X, Y and 
Z. In Ref. 218| . the authors showed that arbitrary rotations of the "spatial 
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parity" , can be manipulated through a series of hnear optical components. 
This allows for the use of the spa tial parity of single photons as a quantum 
bit in quantum information tasks 235| . -For example, Z can be implemented 



by a mirror reflection or a lens system. One can discriminate between even 
and odd spatial modes of a single photon by interfering the field with its mir- 
ror image 236| . This can be achieved using the asymmetric Mach-Zehnder 
interferometer shown in figure [321 Photons traveling through arm 1 of the 



interferometer are subject to 2 mirror reflections, while photons in arm 2 
are subject to 3 mirror reflections. Each reflection transforms the transverse 
wave vector as {qx,<ly) — > {Qx, —%)■ Assuming a single-photon detection at 
an exit port, the photon annihilation operators at the exit ports are related 
to the input port by 



ai(q) 



'Mix, Qy) - r^e*^"ao(gx, -Qy) 



(166) 

aa^qj = itre'^^aoiqx, -Qy) + itre'^^ao{qx, Qy) (167) 

where t and r are the transmissivity and reflectivity of the beam splitters. 
6i and 62 are phases due to propagation through arms 1 and 2. Assuming 
\6i — ^2! mod 27r = 0, and t = r, an even input mode will leave the interfer- 
ometer in output 2 while an odd mode will exit in output 1. The asymmetric 
interferometer is thus capable of projecting onto even or odd modes, and is 



the spatial parity analog to the polarizing beam splitter [218 



To perform rotations in parity space, it is necessary to imprint a variable 



phase on half of the light beam. In Ref. |218l ] it was shown that this can be 
done using two phase plates. Each plate is placed in half of the input field, 
so that plate 1 covers the region x < and imprints a phase exp(— 2^/2), and 
plate 2 covers the region x > and imprints a phase exp(2^/2). The parity 
rotation depends upon the overall phase difference between the two halves of 
the field. 

Using a set of variable phase plates and asymmetric interferometers for 
each of the down-converted fields, Yarnall et al. violated Bell's In equa lity in 



spatial parity space for pump beams with different spatial parity [219 



9. Applications to quantum information 

The facility with which one can produce highly-entangled states from 
SPDC and manipulate the spatial characteristics of the pump and down- 
converted fields allows for the investigation of a wide range of quantum in- 
formation applications and phenomena. A principal advantage of the spatial 
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degrees of freedom of photons, as compared to the polarization degree of 
freedom, for example, is the high- dimensionality of the system. This al- 
lows for the generation of entangled continuous variable states or entangled 
(i-dimensional qudits. Higher dimensional {D > 2) quantum systems have 
been shown to allow fo r hig her security and transmission rates in quantum 
cryptography 2371 . l238l. l239j , increased s ecurity iri two-party protocols such 



as bit commitment 124011 and coin tossing j24lL |242| |. and more robust tests of 



Bell's inequalities 243l |. In sections [7] and [8] we showed that the two-photon 
quantum state produced by SPDC can be highly entangled. Furthermore, 
there exist schemes for engineering the two -photon enta nglement by properly 
shaping the two-photon angular spectrum 16|, l244l . |245| | . There is also a wide 
array of gadgetry devise d to manipulat e these spatial degree s of freed om, in- 
cluding lens systems 21 . 168 . 167 . 246 1. holographic masks ( l84l 247 1, phase 
plates ^209l . l210l . l219| and more recently the development of spatial light mod- 
ulators have opened up even more possibiliti es for the synthesis and control 



of the spatial degrees of freedom of photons |248l . |249| . |250| . |251 



Orbital angular momentum is a natural choice for encoding single- and en- 
tangled qudits, and a number of proof of principle experiments have been per- 
formed which employ th ese d egrees of freedom. These inc lude the demonstra- 
tion of bit commitment 252| and coin tossing protocols 176| with entangled 



qutrits {D = 3), manipulation of heralded qutrits [l75|, violatio n of Bell' s 



173, 174 



inequality and entanglement concentration of entangled qutrits 
violation of Bell's inequality in vari ous two-dimensional subspa ces of th e 
infinite-dimensional 0AM space and quantum cryptography |253l . l254i . 

Higher-dimensional quantum information protocols can also be imple- 
mented using the linear transverse position and momentum. In this respect, 
entangled qudits have been produced using Z)-slit apertures 2^, 255], and by 
binning the near and far-field detection positions 256l | and quantum cryp- 



tography with a 3 7- dimensional alphabet has been demonstrated 257 



Another notable use of spatial entanglement is the creation of mul tiply 
entangled states, which are entangled in multiple degrees of freedom 2581 . 



259!]. In fact, pairs of photons which are hyperentangled (entangled in a ll 



degrees of freedom) have been p roduced and t omographically analyzed 260l 



It h as been shown theoretically [2611 . |262| . 12631 ] and experimentally [ 2011 . 12641 . 
2651 ]. that entanglement in an auxiliary degree of freedom can improve Bell- 
state analysis of photon pairs. Barreiro, Wei and Kwiat have used auxiliary 
0AM entanglement to surpass the channel cap acity imposed by linear optics 
in the quantum superdense coding protocol 2651 ]. First-order transverse 
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La 



Figure 33: Multiple-slit scheme used to produce entangled qudits. 



modes and polarization can be used to define single-pho ton logica l qubits 



which are invariant under azimuthal coordinate rotations 266, 267 , and to 



implement small-scale quantum algorithms [268 



9.1. Generation of Entangled Qudits 

There are a number of ways to generate entangled qudits using the spatial 
degrees of freedom of twin photons. A natural choice is using a set of trans- 
verse modes with discrete quantum numbers, such as the Hermite- Gaussian 
or Laguerre-Gau ssian modes. A number of experiments in this direction have 



been performed |l73L 1 174 1 1751 . 1 1761 . 1252 



Another possibility is through D-slit apertures placed in each down- 



converted field [2551 . |20|. as illustrated in Fig. |33l By focusing the pump 
beam in the plane of the apertures, one can guarantee that the photons al- 
ways pass through symmetrically opposite slits of their respective apertures. 
In this fashion, the two-photon state immediately after the apertures can be 
written as 



20, 255 



(D-l)/2 

J2 exp{ikd^f/2z)\j),\-j),, (168) 

i=(l-D)/2 

where d is the distance between the slits, z is the distance from crystal to 
slits, k is the wave number and \j) are states describing the spatial profile of 
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a single photon which has passed through sht j, which are given by 



\j)= / dpR{p-jd)\p). 



(169) 



Here R{p) are the aperture functions describing each individual slit. The 
position correlation of these photons as well as the entanglement can be 
confirmed by measuring in the near-field and far-field of the apertures, re- 
spectively. In reference 25l| . it was demonstrated that different spatial qu- 



dit states can be generated by sending one of the down-converted beams 
through diffractive apertures with controllable transmission coefficients, pro- 
duced with a spatial light modulator. 

Another technique for generating spatial qudits is by simply defining a 
discret e arr ay of spatial regions. This has been realized in the single-photon 
257l | for dimension D = 37 and entangled-photon regime for up to D = 



re gime 



6 |256|. Using lens systems, it is possible to create and measure superposition 



states 269 



Possible applications of two entangled qudits in quantum communication 
rely on the ability of transmitting these photonic high dimension entangled 
states between different communicating parties. The propagation of entan- 
gled states of qubits have been subject of recent studies. Optical-fibers links 
were used to send energy-time entangled qubits for recei vers separated by 



more than 10km 270 



220 showed that 



and a test of the Bell inequality 
the two-photon state was still in an entangled state. Free-space distribution 
of polarization entangled qub its through the atmosphere was demonstrated 



for a distance of 144km |27l|. The experimental free-space propagation of 
two entangled four-dimensional qud its o r "ququarts" produced by four-slit 
apertures was reported in reference 272| . Free space propagation of the en- 



tangled photons was performed on a laboratory scale. The spatial correlation 
of the propagated state was observed at the image plane of two lenses and 
high-fidelity with its original form was measured. Far from the image plane a 
conditional interference pattern of the propagated state was measured. The 
presence of the conditional interference pattern shows that the state of the 



propagated ququarts is entangled [20|. Several studies of free space propaga- 



tion of single and biphoton fields with 0AM have been performed, and are 
mentioned in more detail in section 19.21 

An important step in future use of entangled spatial qudits is the de- 
velopment of quantum tomography suitable for this degree of freedom or a 
method for the determination of the density matrix of the compound system. 
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In references 2731 . l274j | a technique for determination of the density operator 
was developed for the case of two spatial qubits, generated with two double- 
slits placed in the path of the photon pair. Measurements were performed at 
the image plane and in the focal plane of a lens for different positions of the 
detector in the transverse plane. This meth od c an be used for the tomog- 
raphy of two spatial qubits in a mixed state 275| . In a second method, the 
interference patterns obtained at the intermediate plane between the image 
and the focal plane are used for reco nstruction of the quantum state of two 
spatial qubits |276| and qutrits |277 . 



9.2. Quantum Cryptography 

An interesting application for spatial qudits is quantum key distribution, 
in which the increased dimensionality presents benefits not only in terms 
of the information transmission rate but also the securi ty o f the prot ocol . 
It is straightforwar d to generalize the well-known BB84 |278 || or BBM 279 



protocols to qudits |237l . l238l . l239l |. In this case, it is possible in principle to 
send on average log2 D bits per sifted qudit, while an eavesdropper employing 
an intercept-resend strategy would induce a qudit error rate of Ed = 
which saturates at 1/2 for large D. The increased security is due to the 
increased disturbance, since half the time the eavesdropper measures in the 
wrong basis, and consequently sends the wrong state with a probability of 
{D-l)/D. 

One can achieve the advantages mentioned above using discrete sets of 
transverse modes or employing the continuous position and momentum corre- 
lations present in the down-converted photons. For discrete modes, the quan- 
tum cryptography protocol and theoretical security anal ysis should be nearly 
identical to the general cases of qudits covered in Refs. [23 71 . l238j . Theoret- 



ical proposals and experimental investigations in this directio n have focused 
primarily on the orbital angular momentum degree of freedom 28Cll . 253 , 254 . 

The conjugate position and momentum correlations provide an interest- 
ing scenario from a practical point of view and in addition to the advantages 
mentioned above, have been shown to present particular security features. 
Figure shows an illustration of the main idea. Alice and Bob receive spa- 
tially entangled photon pairs. Independently, they choose randomly between 
two conjugate bases: the near field variable p and far-field variable q. The 
random selection can be performed passively using a 50:50 beam splitter. 
Using optical lens systems (represented by lenses NF and FF), they map 
the near and far-field distributions onto their respective detection planes. 



91 




Figure 34: Quantum key distribution using transverse position and momentum of entan- 
gled photons. Photons are sent to Ahce and Bob's input planes via optical transmission 
systems. The random basis selection can be performed using beam splitters (BS). 



Using an array of single-photon detectors, they register the detection posi- 
tion, and thus infer the variable p or q. Several experimental investiga tions 



of this scheme have bee n performed in the case of entangled photons 281 
and single photons 257 



A first motivation for the study of spatial QKD is the increased infor- 
mation transmission rate per photon pair. In the case that Alice and Bob 
measure in the p basis, the secret key rate A/p is gi ven by th e difference 



between Alice-Bob and Alice- Eve mutual information: 282. 283 



where the mutual information is given by 

I{xi, cca) = H{xi) - H (ccilaJa) , 

and 

H(xi\x2) = — dxidx2P{xi,X2)\nP{xi\x2] 



(170) 

(171) 

(172) 
(173) 
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are the e ntropy an d conditional entropy respectively, for a continuous dis- 



Ltrop ; 
23, 



tribution |213l . I284j |. P{x) is the probability distribution for x, P{xi,X2) is 
the joint probability distribution for Xi,X2 and P(a;i|a;2) is the conditional 
probability distribution for Xi conditioned on measurement of 0:2 • A simi- 
lar expression can be written for the secret key rate Alq corresponding to q 
measurements. Here E stands for the Eve 's rn easurement, and we are con- 
sidering the case of forward reconciliation [282|. Using the above equations, 
the key rate can be put in the form 



AI, = H{p^\p^)-H{p^\E), 



(174) 



and similarly for Alg. For indivi dual attacks, one can consider that Alice, 
Bob and Eve share a pure state |282| . Thus, upon measurement. Bob and 
Eve prepare Alice's phot on in some state, which must satisfy the entropic 
uncertainty relation |285|: 



H{p^\E) + H{q^\qs)>\n7re. 
To establish a secret key requires A/ > 0. This leads to (282'] 
A/ > In Tie - H{q^\qs) - H{p^\pj,) > 0. 



(175) 



(176) 



For any d istribution with variance cr^, the entrop y sat isfies H < In ^l-neo"^ 



2131 . l284j . Thus, the key rate can be bounded by |283 



AI>-\^2A{qMB)^{PA\PB). 



(177) 



where A{x\y) is the variance of the distribution of x given a measurement 
of y. 

This condition is fulfilled when A(q^|qr^)A(p^|p^) < 1/2, which is equiv- 
alent to the EPR variance criterion (11121) . For photons produced by SPDC, 
the EPR criteria depends on the characteristics of the pump beam and the 
length of the non-linear crystal L. 

Implementation of QKD requires preservation of quantum information 
over long transmission distances, and as such is a technical challenge. QKD 
with spatial variables should be most aptly suited to a free-space transmis- 
sion, as opposed to optical fibers, due to the multi-mode character of the 
down-converted fields. There are at least two main factors which limit the 
distance at which spatial QKD could be a viable endeavor. The first are fluc- 
tuations in the index of refraction along the free-space transmission path. 
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This causes local phase fluctuations which distort the single-photon fields, 
increasing the error rate. Atmospheric turbulence causes fiuctuations in the 
index of refraction which transform an initial field ip{p) — > exp[z5(p)]'?/^(p), 
where the position-dependent phase 5{p) describes the phase fiuctuations in 
the wave function. The transmission of spatial q uantum in formation un- 
der these conditions has been considered in Refs. 
Laguerre-Gaussian modes and in Ref. 



286, 287 in the case of 



269l | for single-photon spatial QKD. 



Assuming that the phase fiuctuations are isotropic and a random Gaussian 
process, the induced error rate is dependent upon the initial fiel d and th e 



length scale of the phase fiuctuations known as the Frie d par ameter |288l . |286 



Typical values of the Fried parameter are 10 — 60 cm |289| . For realistic ex 



perimental parameters, it was shown in Ref. j269j that spatial QKD should 
be feasible up to a few kilometers. It is important to note that the recent 
progress in the field of classical optical communication has shown that it 
should be possible to correct the effects of turb ulence in real time with adap- 
tive optics and a reference beam 289l |290| . l29lj . A recent experiment by Pors 



et al. has shown that 0AM entangle ment should resist typical atmospheric 



turbulence for up to a few kilometers [292 



The second limiting factor is the presence of losses due to absorption, 
which reduces the transmission rate as well as the signal to noise ratio, which 
Eve can exploit to her advantage. The noise present in spatial QKD comes 
from s purio us ambient photons and detector dark counts. Following Zhang 



et al. |283| . three types of events can contribute to the coincidence count 
rate, occurring with probabilities P(spdc,spdc), P(noise, spdc) or P(spdc, 
noise), and P(noise, noise). Focusing only on dark counts, they show that 
for reasonable experimental parameters, one can tolerate losses of a few dB 
and still satisfy Eq. f ll77p . This level of loss corresponds to free-space trans- 
mission distances on the order of a few km. However, it is important to note 
that the left-side of inequality fll77p is an upper-bound which is saturated 
only when Alice and Bob's coincidence distributions are Gaussian. However, 
assuming a constant dark count rate for each detector, the coincidence counts 
in the presence of loss will be composed of a Gaussian distribution of cor- 
related photons, combined with a constant background count, resulting in a 
non-Gaussian distribution. The interesting eavesdropping strategy for Eve 
is to perform measurements so that she can mask her disturbance with the 
counts due to noise, so that, in this case, her disturbance should appear as 
a constant background noise. This can be achieved with a simple intercept 
resend strategy, so that when she measures in the same basis as Alice and 
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Figure 35: Illustration of teleportation of the spatial degrees of freedom of a single pho- 
ton. A nonlinear crystal is used to perform second harmonic generation(SHG). Bob's 
momentum shift can be performed using spatial light modulation (SLM). 



Bob, she goes undetected. Zhang et al. 283| have shown that consider- 
ing the non-Gaussian nature of the total coincidence distribution allows one 
to increase the allowable losses by at least a factor of ten, giving a secure 
transmission distance of several tens of kilometers. Thus, in a realistic sce- 
nario, the noise due to dark counts can increase the allowable losses, while 
maintaining security of the protocol. 

9.3. Quantum Teleportation 

Quantum teleportation is a method to transfer quantum information be- 
tween distant locations, without direct transmission of the physical system 
containing the information. The original quantum teleportation proposal 
293| considered discrete D-\eYe \{D = 2) syste ms. Further proposals general- 



ized it to continuous variables |294j . l295l . l296j . Quantum teleportation plays 



an important role in quant um information processing since it is the basis 
for q uantum repeaters 297l |. non-local interactions between light and matter 



298| . as well as linear optics quantum computing |299 

A quantum teleportation protocol which us es transve rse spatial properties 
of intense light fields has also been proposed 300l . 301 |. and a teleportation 
procedure which can t ransm it the angular spectrum of a single-photon has 
been proposed in Ref. [302 . 



We will describe the teleportation procedure of Ref. 3021 ] in more detail. 
Figure [33 illustrates the main idea. Alice is in possession of a single photon 
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described by the quantum state 



J dqu{q)\q),, (178) 



where u{q) is the normahzed angular spectrum of the field in the z = plane. 
She would like to send the angular spectrum u{q) of the single-photon field 
to Bob. To do this, Alice and Bob share a pair of spatially entangled photons 
described by the quantum state 



\iPab) oc / / dqidq,v{qi + q^hiqi - qj, (179) 




where v is the angular spectrum of the laser and 7 is the phase matching 
function defined in Eq. fISUl) . Assuming that Alice is in possession of fields 1 
and 2 and Bob of field 3, she follows the usual recipe for quantum teleporta- 
tion 2931]. First, it is necessary for Alice to perform a joint measurement on 
her fields, and sends the measurement result to Bob using classical commu- 
nication. Depending on the outcome. Bob performs a unitary operation on 
his field, and then recovers the initial quantum state of field 1. 

Alice can perform the joint measurement on the spatial degrees of freedom 
of photons 1 and 2 by injecting then into a nonlinear up-conversion crystal cut 
for second harmonic generation (SHG) or sum frequency generation. T he ful l 



multimode Hamiltonian operator for SHG has been described in Ref. [303 
and in the thin-crystal limit can be written as 



Hoc dqdq'dq"6\q + q' -q")aiiq)a2iq')aliq") + H.c.. (180) 




Index 4 denotes the second-harmonic field, which is assumed to be initially 
in the vacuum state. After the nonlinear interaction, Alice measures the 
transverse wave vector q" of the up-converted field, and sends the measure- 
ment result to Bob. The detection position is sent to Bob by a classical 
communication channel, as shown in figure [35l The paraxial field operator 
which describes photodetection at position in the far-field is given by Eq. 
(11081) . Applying this operator to the state at the output of the beam splitter, 
detection of a photon at the position p^, selects the state 

|V^/)(X 11 dqdq'u{q)<^(^jPo-q,q'^ \q'),, (181) 
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where all modes in the vacuum state have been excluded, n is the wave 
number of the SHG field, / is the focal length of the lens, and $ is the 
two-photon angular spectrum defined in Eq. fl57|) . 

In the limiting case of an EPR-like state, in which the angular spectrum of 
the SPDC pump field is a plane wave, v{q) oc 6{q), and the phase matching 
function is very broad, so that it is reasonable to approximate 7(9) oc 1. 
This is the thin crystal approximation described in section 13.41 Under these 
conditions the state teleported to Bob fllSip is 



\^r") = j dqu[q+-^pj,^\q),, (182) 

which is equivalent to the original state ( I178p except for a translation in q. 
Applying a momentum translation: 

P(a)|g) = |g + a), (183) 

Bob will then have P{Kpj^/ f) \'^f^^) = 10)3, which corresponds to perfect 
teleportation of the initial state \4>)i- The momentum translation can be 
implemented using a spatial light modulator. 

In a more general case where 7 cannot be approximated by unity, the 
final state after momentum translation is 

\^f) =C jj dqdq'u{q)viq' - 9)7 (^2jPo - q' - q^ \q), , (184) 

where C is a normalization constant. Assuming that Alice detects a photon 
in b eam 4, the fidelity measures the success of the teleportation protocol 



202| . and is defined as J-" = \{(f)\^ljf)\. The fidelity depends on the format of 
the angular spectrum v{q) and phase matching function 7(9), respectively. 

This teleportation procedure has not yet been realized experimentally, 
and prevents several serious technical challenges. The principal obstacle is 
the low up-conversion efficienc y for single-photon fields. New materials, such 
as periodically-poled crystals [304] or photonic crystal structures may im- 
prove these conversion rates. In spite of the low conversion efficiency, a 
similar experiment using pairs of non-linear crystals to implement a com- 
plete Bell-state measurement in the polarization degree of freedom has been 



performed [305 



The same scheme can be used to swap entanglement between spatially 
entangled photons pairs, which would be necessary for long-distance spatial 
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quantum communication using a large portion of the available Hilbert space. 
Another application of this teleportation procedure is the frequency tuning 
of one photon of a spatially entangled pair. This can be done by replacing 
photon 1 with an intense field, increasing the up-conversion rate. Frequency 
tun ing o f entangled photons has been realized using a periodically-poled crys- 



tal 304 



10. Conclusion 

Parametric down- conversion has definitely played a very important role 
in the development of Quantum Optics and more recently Quantum Infor- 
mation. This review focused on the spatial degrees of freedom of the down- 
converted fields, which have been an important ingredient. The quantum 
correlations between down-converted photons have allowed for the produc- 
tion of quantum states of light, leading to the observation of non-classical 
effects and the study of entangled states and their application in Quantum 
Information tasks. Initially, the concept of spatial coherence in Optics was 
generalized to include two-photon coherence. Counter-intuitive aspects of 
the quantum world have been experimentally observed with the transverse 
spatial degrees of freedom, such as the measurement of the De Broglie wave- 
length of a two-photon wave packet, the anti-bunching of photon pairs, condi- 
tional two-photon interference patterns, quantum erasure, ghost images and 
ghost interference. These investigations have led to a considerable concep- 
tual improvement and a real advance in the knowledge of the quantum world 
through the study of the spatial correlations in parametric down-conversion. 
Applications immediately arose from these new concepts, such as quantum 
lithography and metrology, quantum cryptography schemes, quantum tele- 
portation and spatial qudits for quantum computation. 

This field is still very active and contributions in the fundamental and ap- 
plied domains appear regularly in the literature. Among the most promising 
research in the field, we would like to mention the study of the transverse cor- 
relations as a resource for quantum information in the domain of continuous 
variables and high dimension discrete states. Transverse spatial variables are 
natural entangled continuous variables that can be prepared in a vast range 
of entangled states whose properties and applications are just beginning to 
be explored. 
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